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Abstract We review the theory of relaxational kinetics of superfluid turbulence — 
a tangle of quantized vortex lines — in the limit of very low temperatures when the 
motion of vortices is conservative. While certain important aspects of the decay 
kinetics depend on whether the tangle is non-structured, like the one correspond- 
ing to the Kibble-Zurek picture, or essentially polarized, like the one that emu- 
lates the Richardson-Kolmogorov regime of classical turbulence, there are com- 
mon fundamental features. In both cases, there exists an asymptotic range in the 
wavenumber space where the energy flux is supported by the cascade of Kelvin 
waves (kelvons) — ^precessing distortions propagating along the vortex filaments. 

At large enough wavenumbers, the Kelvin-wave cascade is supported by three- 
kelvon elastic scattering. At zero temperature, the dissipative cutoff of the Kelvin- 
wave cascade is due to the emission of phonons, in which an elementary process 
converts two kelvons with almost opposite momenta into one bulk phonon. 

Along with the standard set of conservation laws, a crucial role in the theory of 
low-temperature vortex dynamics is played by the fact of integrability of the local 
induction approximation (LIA) controlled by the parameter A = ln(A/ao), with 
X the characteristic kelvon wavelength and ao the vortex core radius. While ex- 
cluding a straightforward onset of the pure three-kelvon cascade, the integrability 
of LIA does not plug the cascade because of the natural availability of the kinetic 
channels associated with vortex line reconnections. 

We argue that the crossover from Richardson-Kolmogorov to the Kelvin-wave 
cascade is due to eventual dominance of local induction of a single line over the 
collective induction of polarized eddies, which causes the breakdown of classical- 
fluid regime and gives rise to a reconnection-driven inertial range. 
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1 Introduction 

Superfluid turbulence (ST), also known as quantum/quantized turbulence, is a tan- 
gle of quantized vortex lines in a superfluid lIL 12, |j, |4i] . In the last one and a 
half decade, and especially in the recent years, the problem of zerojioint ST has 
evolved into a really hot subfield of low-temperature physics lialdlTLIsLlgLfl^ fllL 
El HHIil, [1113, [Sin, 113, Hi, 22, 23, mm HI. In this work, we summa- 
rize theoretical developments of the two authors on the theory of decay of ST at 
T = 0, previously published in a number of short papers/letters 11131 [l4l[l5ll2ll 12311 . 
To render the discussion self-contained, we also review the analysis of Ref. ||5|], 
where a number of conceptually important facts about zero-point ST has been 
revealed. 

Superfluid turbulence (ST) can be created in a number of ways: (i) in the 
counter-flow of normal and superfluid components 1271 l28h (ii) by vibrating ob- 
jects 10, [3 [13, [H], (iii) as a result of macroscopic motion of a superfluid (re- 
ferred to as quasi-classical turbulence), in which case ST can mimic, at large 
enough length scales, classical-fluid turbulence |[l^ H9lll[30l[3ll[T8l[l9llp|, Ov) 
in the process of (strongly) non-equilibrium Bose-Einstein condensation 11321 l33ll . 
in which case it is a manifestation of generic Kibble-Zurek effect ll34ll . 

During about four decades — from mid fifties till mid nineties — ST turbulence 
was intensively studied in the context of counter-flow of normal and superfluid 
components. An impressive success has been achieved in the field, with the promi- 
nent contributions by Vinen — the equation for qualitative description of growth/ 
decay kinetics (Vinen's equation), extensive experimental studies ll27ll . — and pi- 
oneering microscopic simulations of vortex line dynamics by Schwartz fz^. The 
counter-flow setup naturally implies finite density of the normal component, and 
thus a relatively simple — by comparison with the T = case — relaxation mecha- 
nism. The normal component exerts a drag force on a vortex filament. As a result, 
the vortex line length decreases. For an illustration, consider a vortex circle of 
the radius R. In the absence of normal component, the ring moves with a con- 
stant velocity V = V{R), the radius R remaining constant. With the drag force, the 
ring collapses, obeying a simple law R oc —aV{R), where a is the dimensionless 
friction coefficient measuring the drag force in the units of Magnus force. Sim- 
ilarly, the drag causes the decay of Kelvin waves — precessing distortions on the 
vortex filaments — and the value of a~' gives the number of revolutions a dis- 
tortion makes before its amplitude significantly decreases. At a ^ 1, the decay 
scenario of a (non-structured) vortex tangle is as follows. The vortex lines recon- 
nect producing the distortions (Kelvin waves), these distortions decay due to the 
drag force thereby reducing the total line length and rendering the tangle more and 
more dilute. 

In contrast to the counter-flow setup, where the case a <C 1 can be considered 
rather specific, for the non-equilibrium process of Bose-Einstein condensation in 
a weakly interacting gas, the regime a ^ 1 is quite characteristic in view of the 
following circumstance. The condensation kinetics is a classical-field process the 
essential part of which is quantitatively captured by the time dependent Gross- 
Pitaevskii equation (a.k.a. non-linear Schrodinger equation). Hence, the most nat- 
ural setup is when the initial condition for the process features large occupation 
numbers for the bosons, and the whole process can be accurately described by the 
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Gross-Pitaevskii equation from the very outset (see Fig. [T] for an illustration). In 
the idealized situation of the classical field, the final temperature asymptotically 
approaches zero in view of the ultraviolet catastrophe. Another situation where 
the regime a <C 1 occurs naturally is the quasi-classical turbulence, which, unlike 
counterflow turbulence, is due to a macroscopic turbulent motion and is generated 
by purely classical means at arbitrary temperatures without the help of mutual 
friction. 

At a <C 1 (a °= at r ^ 0, see Ref. (35']), dynamics of vortex lines in a 
significant range of length scales become essentially conservative. In this case, re- 
laxation of turbulent motion in non-linear systems usually involves certain types 
of cascades. In both classical and superfluid turbulence the key role is played by 
the cascade of energy in the wavenumber space towards high wavenumbers. A 
realization of this peculiar relaxation regime is due to the following conditions 
satisfied by the system: (i) Aen ^ ■^cutoff — a substantial separation of the energy- 
containing scale Aen from the scale Acutoff where the dissipation becomes appre- 
ciable, which defines the cascade inertial range Aen > A > Acutoff- (ii) the kinetics 
are local in the wavenumber space, i.e. the energy exchange is mainly between 
adjacent length scales, and (iii) the "collisional" kinetic time Tcoii(A) — the time 
between elementary events of energy exchange at a certain scale A — gets progres- 
sively shorter down the scales. These requirements determine the main qualitative 
features of the cascade: The decay is governed by the slowest kinetics at Aen, 
where the energy flux £ in the wavenumber space is formed, while the faster ki- 
netic processes at shorter scales are able to instantly adjust to this flux supporting 
the transfer of energy towards Acutoff, where it is dissipated into heat. Thus, the 
cascade is a (quasi-)steady-state regime in which the energy flux e is constant 



4 



through the length scales and the variation of £ in time happens on the longest 
time scale ^ Tcoii(/len)- 

With vortex lines, one can think of quite a number of different cascades. Per- 
haps the most obvious one is the Richardson-Kolmogorov cascade of eddies, 
which can be realized in a superfluid (even at T = 0) due to its ability to emu- 
late the classical-fluid turbulent motion by the motion of a polarized vortex-line 
tangle. Clearly, this type of the cascade is fundamentally impossible in the non- 
structured Kibble -Zurek-type tangle. For quite a long time it was believed that 
non-structured ST decays at T = through the Feynman's cascade of vortex rings 
13^ . Feynman conjectured that reconnections of the vortex lines produce vortex 
rings, with subsequent decay of each ring into a pair of smaller rings, and so forth. 
It can be shown, however, that this conjecture is inconsistent with simultaneous 
conservation of energy and momentum (see Ref. Isj] and Sec.|4ll. Another obvious 
option is a cascade of Kelvin waves supported by non-linear interaction (scatter- 
ing) of kelvons. While such a cascade is indeed possible (see Ref. fl3h and Sec.[3]l, 
it is subject to a peculiar constraint on the maximal energy flux. The origin of this 
constraint is the approximate integrability of the vortex dynamics controlled by 
the large parameter 

A=ln(A/flo), (1) 

where A is the characteristic kelvon wavelength and ciq is the vortex core radius. 
If the parameter A is large, then the leading term of the Kelvin-wave dynamics is 
given by the so-called local induction approximation (LIA) in which the velocity 
of an element of the vortex line is due to local differential properties of the line. 
The LIA dynamics turns out to be integrable, so that mere non-linearity of the LIA 
equations of motion does not lead to the Kelvin-wave cascade. This circumstance 
predetermines the existence of intermediate cascades that are necessary to transfer 
the energy to the wavenumbers high enough for the pure Kelvin-wave cascade to 
take over. All the intermediate cascades have to involve reconnections, to lift the 
constraint imposed by the integrability of LIA. In the non-structured tangle, there 
is only one type of reconnection-induced cascade. An elementary step of this cas- 
cade at the wavelength A consists of two adjacent events: (i) emission of a vortex 
ring with the radius ^ A by a local self-crossing of a kinky vortex line and (ii) 
re-absorbtion of the ring by the tangle. Both events are accompanied by transfer- 
ring Kelvin-wave energy to shorter wavelengths (but still on the order of A), since 
reconnections of the vortex loops produce Kelvin-wave structure with wavelength 
smaller than the loop radii. For the self-reconnection cascade to be efficient, the 
vortex line has to be kinky (loosely speaking, fractalized) at all relevant length 
scales. Due to the integrability of LIA, the fractalization is provided by the cas- 
cade itself: In the absence of self-crossings at a given length scale, the amplitude 
of Kelvin waves keeps growing due to self-crossings at larger length scales. It is 
worth emphasizing that in contrast to Feynman's scenario, where the vortex rings 
are the energy carriers, in the local self-crossing scenario, the rings play only a 
supportive role. In the process (i) a ring is just a product of self-crossing and in 
the process (ii) the ring is just a cause of yet another reconnection event. Clearly, 
the crucial part is played by Kelvin waves, which carry the energy, and reconnec- 
tions, which directly promote the Kelvin-wave cascade. 

In the polarized tangle emulating the Richardson-Kolmogorov classical-fluid 
cascade, there are two more varieties of reconnection-supported Kelvin-wave cas- 
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cades. One is when reconnections are between the bundles of quasi-parallel vortex 
lines, and the other one when the reconnections are between two neighboring lines 
in a bundle. In Sec.[5]we present the arguments ll2lll that in the theoretical limit of 
A ^ oo, all the three reconnection driven cascades are necessary to cross over from 
Richardson-Kolmogorov to pure Kelvin-wave regime, the crossover being a series 
of three distinct cascades: bundle driven, neighboring reconnection driven, and lo- 
cal self-crossing driven. The total extent of the crossover regime in the Kelvin 
wavenumber space is predicted to be A/10 decades. With realistic A < 15 (for 
^He), this means that the crossover takes about one decade, within which one can 
hardly expect sharp distinctions between the three different regimes. 

At r = 0, the reconnection-supported cascade(s) ultimately crosses over to the 
pure Kelvin-wave cascade, and the latter is cut off by phonon emission Isl flOlflSl] . 
The elementary process of phonon emission converts two kelvons with almost 
opposite momenta into one bulk phonon ifish . 

The rest of the paper is organized as follows. In Sec.|2]we render basic equa- 
tions and notions of vortex dynamics: Biot-Savart equation and the LIA, Hasimoto 
representation for the LIA in terms of curvature and torsion, conservation of en- 
ergy and momentum, Hamiltonian formalism for the Kelvin waves, conservation 
of angular momentum, implying conservation of the total number of kelvons. We 
pay special attention to the problem of ultraviolet regularization of the theory, 
crucial for an adequate treatment of non-local corrections to the LIA. In Sec. |3] 
we develop a theory of the pure Kelvin-wave cascade (supported by elastic scat- 
tering of three kelvons — the leading non-trivial scattering event consistent with 
conservation of energy, momentum, and the total number of kelvons). We show 
that the leading term in the three-kelvon scattering amplitude does not contain lo- 
cal contributions, in accordance with the integrability of LIA. We then find the 
spectrum and energy flux of the pure cascade. In Sec.|4]we analyze the scenario 
of local self-crossings, starting with the explanation why pure Feynman's cascade 
cannot be realized in view of conservation of momentum and energy. We show 
that this scenario implies fratalization of the vortex lines, and derive correspond- 
ing spectrum of Kelvin waves. We then discuss a simple Hamiltonian model jsj] 
featuring a collapse-driven cascade in otherwise integrable system, and argue why 
the spectrum of the cascade in this model is similar to the one in the scenario 
of local self-crossings. In Sec.|5]we address the problem of the crossover from 
Richardson-Kolmogorov to Kelvin-wave cascade. Sec.|6]is devoted to the theory 
of kelvon-phonon interaction. Starting from the hydrodynamical Lagrangian, we 
derive the Hamiltonian of kelvon-phonon interaction. The Hamiltonian allows us 
to straightforwardly formulate kinetics of the kelvon-phonon processes, and, in 
particular, find the cutoff momentum for the pure Kelvin-wave cascade. In Sec.|7] 
we summarize the main qualitative aspects of the theory and put our findings in 
the context of experiment and simulations, including the ones that are still miss- 
ing. We also formulate certain theoretical questions to be addressed in the fu- 
ture. Finally, we critically discuss the concept of bottleneck between Richardson- 
Kolmogorov and Kelvin-wave cascades, put forward recently in Ref. ll20h (see also 
Ref. 122]), which we disagree with on the basis of our theoretical analysis. 
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2 Basic relations 

2.1 Biot-Savart law. Energy and momentum 

The zero-temperature hydrodynamics of a superfluid is nothing but the hydrody- 
namics of a classical ideal fluid with quantized vorticity — the only possible rota- 
tional motion is in the form of vortex filaments with a fixed value of velocity cir- 
culation K. Therefore, one can apply the Kelvin-Helmholtz theorem of classical 
hydrodynamics — stating that the vortices move with the local fluid velicity — to 
obtain a closed dynamic equation for vortex filaments in a superfluid. Alterna- 
tively, one can start with the hydrodynamic action for the complex-valued field 
and derive the equation of vortex motion from the least-action principle. In the 
present section, we use the former approach, readily yielding the answer. [The lat- 
ter approach will be used in Sec. 16. 11 where it will prove crucial for describing the 
vortex-phonon interaction.] 

If the vortex lines are the only degrees of freedom excited in the fluid, and if the 
typical curvature radius and interline separations are much larger than the vortex 
core size, ao> then the instant velocity field at distances much larger than aq from 
the vortex lines is defined by the form of the vortex line configuration. Indeed, 
away from the vortex core the density is practically constant and the hydrodynamic 
continuity equation reduces to 

V • V = . (2) 

Then, taking into account that V x v = everywhere except for the vortex lines, 
and that for any (positively oriented) contour F surrounding one vortex line 

\-dl=K, (3) 

r 

we note a direct analogy of our problem with the magnetostatic problem of finding 
magnetic field produced by thin wires: Velocity field is identified with the mag- 
netic field, and the absolute value of the current of each wire is one and the same 
and is proportional to K. The result is given by the Biot-Savart formula 

^' 4nJ |r-sP ' 

where radius- vector s runs along all the vortex filaments. 

According to Kelvin-Helmholtz theorem, each piece of the vortex line should 
move with a velocity corresponding to the velocity of the net motion of small 
contour (of the size, say, of order ciq) surrounding this element. This fact is very 
important, since the velocity field ^ is singular at all points on the vortex line, and 
the Kelvin-Helmholtz theorem yields a simple regularization prescription: Take a 
small circular contour with the center of the circle at some vortex line point and the 
plane of the circle perpendicular to the vortex line (at the point of intersection), and 
average the velocity field around the contour — to eliminate rotational component 
of the contour motion. 

As is clear from the singularity of the velocity field at some point s = Sq in 
the vortex line comes from the integration over the close vicinity of the point sq. 
To isolate the singularity, we expand the function ) around the point sq = s(^o): 

s = so + s'^^+s'^^^V2 + -.- • (5) 
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Here ^ is the parameter of the line. It is convenient to use the natural parameteri- 
zation, that is to choose E, to be the (algebraic) arc length measured from the point 
Sq. Substituting this expansion into (|4|, we get 

v(r->so) = g^s'^xs'^'^ j ^ —+ regular part, (6) 

where is some upper cutoff parameter on the order of the curvature radius at the 
point So- The integral in the right-hand side is divergent at ^ ^ 0. To regularize 
it we note (i) that, by its very origin, the expression ^ is meaningful only at 
I I > ao, and (ii) the part of the vortex line with | ^ | <C ao does not give a significant 
contribution to the net velocity on the contour of the radius ^ ao- Hence, with a 
logarithmic accuracy we can adopt the regularization \^\ > ao, that is 

dE 

^ ^ -> 21n(^,/ao) . (7) 

Moreover, by fine-tuning the value of aq in accordance with a model-specific be- 
havior at the distances ^ ao, the logarithmic accuracy of the regularization (|7]l can 
be improved to the accuracy ^ ao/|*- We discuss this option in detail in Sec. l2.4l 
and utilize in Sec.[3j 

We thus arrive at the Biot-Savart equation of vortex line motion 

K f dSQ X (S - So) 

'-4^.1 |s-sop ' 

with the integral regularized in accordance with (|7]l. 
Equation ^ has two constants of motion: 

<^=/;^, (9) 

./ |s-so| 

^ = fsxds, (10) 



which (up to dimensional factors) are nothing than the energy and momentum, 
respectively. 



2.2 Local induction approximation (LIA) 

In the absence of significant enhancement of non-local interactions by polarization 
of the vortex tangle, the regular part in (|6]l is smaller than the first term containing 
large logarithm. In such cases it is often — but not always (!), see the theory of 
the pure Kelvin-wave cascade — safe to neglect the second term in ^ and proceed 
within LIA: 

s = j3s'^xs'^'^, (11) 

where 

l5 = ^ln{R/ao), (12) 
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with the typical curvature radius 7? treated as a constant. Recalling that the param- 
eter ^ in (|6ll is the arc length, it is crucial that equation of motion (1111) is consistent 
with this requirement. [Speaking generally, in the course of evolution might de- 
viate from the arc length.] The consistency is established by directly checking that 
(II lb implies 

4- ^ds • t/s = 0. (13) 

at 

From Eq. ( fT3] l it trivially follows — by integrating over ^ — that the total line length 
is conserved. The total line length in LIA plays the same role as the energy (|9| 
in the genuine Biot-Savart equation. Indeed, from (|9]l it is seen that within the 
logarithmic accuracy S is proportional to j3 times line length. It is not difficult to 
also make sure that Eq. dl lb conserves 0^ dlOb . 

The standard constants of motion — the energy (line length), momentum, and 
angular momentum — are not the only quantities conserved by Eq. dl lb . For exam- 
ple, the integral of the square of the curvature radius is also conserved jjTh : 

In the next section we will see that the constant of motion dl4b is just one of the 
infinite set of constants of motion implied by the integrability of LIA. 



2.3 Hasimoto representation. Integrability of LIA 



Betchov OTIl revealed certain interesting properties of LIA, e.g., Eq. (fT4] |. by re- 
writing Eq. dl lb in terms of intrinsic variables of the vortex line: curvature, i^, and 
torsion, T. Hasimoto |!38] further advanced these ideas by discovering that for the 
complex variable ,/), such that 

^ = Wl ^ = -^ (15) 

{•^ is the phase of i//) the LIA equation (fTTb is equivalent to the non-linear Schrodinger 
equation (time is measured in units j3 ~ ^ ) 

The one-dimensional non-linear Schrodinger equation is known to be an inte- 
grable system featuring and infinite number of additive constants of motion, the 
explicit form of which is given by [39] 



In = j (pn{^)d^, 



<Pi = ^IV^l', (17) 
d 



ni+f?2=n 



The integrability of LIA renders any non-trivial relaxation kinetics impossible, 
unless the reconnections are involved to break the conservation of /,,'s. 
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2.4 Hamiltonian formalism 



Suppose a vortex line can be parametrized by its Cartesian coordinates x and y 
as single-valued functions of the third coordinate z. In this case, the Biot-Savart 
law can be cast into a Hamiltonian form, very convenient for our purposes. First, 
introduce a vector p(z,/) = {x{z,t) ,y{z,t)) An contrast to the earlier-defined vector 
s(^,f) that follows the motion of the element of fluid containing the element of 
vortex line, the vector p (zq, /) just defines the geometrical point of intersection of 
the vortex line with the plane z = Zo- With this distinction in mind, it is easy to 
relate the time derivatives of the two vectors: 

Here z is the unit vector in the z-direction. Substituting for 5s /dt the right-hand 
side of Eq. expressing then s in terms of p, and finally replacing p with a 
complex variable w{z^ t) = x{z,t) + iy{z,t), we arrive at the Hamiltonian equation 
of vortex line motion: 

(19) 

ow* 



[\+RtW*{z\)w'{z2)]dz\dz: 



H = ( k/Ak) / 7, .r'" ■ (20) 



\/{z\ -Z2)^+ |w(zi) - W(Z2)P 

The Hamiltonian ( 120b is singular at zi Z2 and thus needs to be regularized. To 
this end we introduce r* such that 

flo<'-*<A, (21) 

and write 

//=//loc+//n.l. + ^(r*/A), (22) 

\+RQw'*{zi)w'{z2)]dzidz2 



Hn.i = (K/An) I 



-Z2)^ + |h'(zi)-w(z2)P 

\-Z2\>rt 



(23) 



//loc = 2/3 y dz^l + \w'iz)\^, (24) 

j3 = ()c/27r)ln(r,/fl,)- (25) 

Here the value of a* ^ ao is specially tuned to eliminate a factor of order unity in 
the logarithm. — The absence of such a factor in Eq. (125 b should not be confused 
with a lack of control on first sub-logarithmic corrections. In view of the first in- 
equality in (|2T|, the Hamiltonian H^.i. is of purely hydrodynamic nature, while the 
Hamiltonian Hioc takes care of both hydrodynamic and microscopic (and system 
specific) features at distances smaller that r*. In view of the condition A ^ ao, im- 
plied by Eq. (121b . the microscopic specifics of the system is completely absorbed 
by the proper value of . Indeed, whatever is the physics at the distances of the 
order of the vortex core radius, the leading contribution from these distances to the 
energy of a smooth vortex line should be directly proportional to the line length; 
and this is precisely what is expressed by Eq. 
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The freedom of choosing a particular value of within the range (1211) can be 
used to introduce the local induction approximation by requiring that 

ln(r*/ao) >ln(A/r*). (26) 

In this case, the Hamiltonian (122b . with 

l5 = {K/2K)ln{X/ao)=AK/2n, (27) 

captures the leading — as long as the interline interactions are not relevant — logarithmic 
contribution, compared to which the non-local Hamiltonian ( 123b can be omitted 
with logarithmic accuracy guaranteed by the parameter A ~ ^ <^ I. 
Less obvious technical trick is to formally set 

= fl* (28) 

to nullify the Hamiltonian //loc- Doing so might seem to violate the range of appli- 
cability of the essentially hydrodynamic Hamiltonian H^ i fiol]. Nevertheless, it is 
readily seen by inspection that under the condition A ^ aq the resulting theory is 
equivalent to the theory (I22ll-( l25l l up to negligibly small corrections of the order 
cio/Pi. The resulting Hamiltonian is 

J V(Z1 -^2)^ + |w(zi) - W(Z2)P 

\Z]~Z2\>a* 

In direct analogy with the pseudo-potential method in the scattering theory, it has 
a status of a piCMiio-Hamiltonian in the sense that, being a convenient but rather 
inadequate model at the scales ^ ciq, it accurately accounts for the long-wave 
motion of the vortex lines (including all sub-leading corrections to LIA coming 
from the distances ^ ao) by an appropriate choice of a*. In the next subsection, we 
describe an explicit procedure of extracting the value of a* from a model-specific 
dispersion relation of Kelvin waves. 



2.5 Kelvons. Angular momentum and the number of kelvons 

Suppose the amplitude of Kelvin waves (KW) is small enough, so that the follow- 
ing condition is satisfied; 

0C[ZUZ2) = . . <1. (30) 

\Zl -Z2\ 

Linearizing ( 129b up to the leading orders in the dimensionless function (X{z\,Z2), 
we obtain 



Ho = (K/87i) f ^£i^[2Rew'*(ziy(z2)-a']. (31) 

J \Zl -Z2 
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The Hamiltonian Hq describes the linear properties of KW. It is diagonalized by 
the Fourier transformation w{z) = L~^^^Y,kWk^'^~ is the system size, periodic 
boundary conditions are assumed): 

Ho = {K/4n)'£o),wlwk, (32) 

k 

yielding Kelvin's dispersion law 



K 

Oik = :7- 
4n 



k\ (33) 



with 



Co = 2 /" — [cos X + (cosx — 1 ) A^] + 

2 / — [cosx - 1/2 + (cosx - \)/x^] w -2.077. (34) 
Jo X 



Generally speaking, the term (? {(ka-tY) should be omitted, since it goes beyond 
the universality range of the pseudo-Hamiltonian ( 1291 ). — In each particular micro- 
scopic model this term will be sensitive to the details of the physics at the length 
scale ^ ciQ- [With logarithmic accuracy, the constant Cq can be omitted as well, 
while the specific value of can be replaced with just and order-of-magnitude 
estimate ao-] 

The practical utility of Eq. ( 133b with sub-logarithmic correction Co ( 134b for a 
given specific microscopic model — different from the pseudo-Hamiltonian ( |29b — 
is as follows. Eqs. (I33b -(l34b can be used to calibrate the value of a*, and thus fix 
the pseudo-Hamiltonian, by separately solving for the KW spectrum in the given 
microscopic model and then casting the answer in the form (I33b -( l34b . Likewise, 
for a realistic strongly-correlated system like ^He the appropriate value of a* in the 
pseudo-Hamiltonian (129 b can be calibrated by an experimentally measured kelvon 
dispersion. 

Although the problem of KW cascade generated by decaying superfluid turbu- 
lence is purely classical, it will be convenient to approach it quantum mechanically- 
by introducing KW quanta, kelvons. In accordance with the canonical quantiza- 
tion procedure, we understand Wk as the annihilation operator of the kelvon with 
momentum k and correspondingly treat w(z) as a quantum field. The Hamiltonian 
functional ( 129b is proportional to the energy — with the coefficient Kp /2, where 
p is the mass density — but not equal to it. This means that if one prefers to work 
with genuine Quantum Mechanics rather than a fake one (for our purposes, the 
latter is also sufficient), using true kelvon annihilation operators, dj^, field operator 
w, and Hamiltonian, H, he should take into account proper dimensional coeffi- 
cients: H = {Kp/2)H, Wk = ^2h/Kpdk. By choosing the units h = K = \, p = 2, 
we ignore these coefficients until the final answers are obtained. 

In the quantum approach, there naturally arises the notion of the number of 
kelvons. This number is conserved in view of its global U{\) symmetry of the 
Hamiltonian, w e'^w reflecting the rotational symmetry of the problem. The 
rotational symmetry is also responsible for the conservation of the angular mo- 
mentum component along the vortex line direction. Thus the number of kelvons is 
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immediately related to the angular momentum component ll4]]l : one kelvon carries 
a single (negative) quantum, —Ti, of the angular momentum component along the 
vortex line relative to the macroscopic angular momentum of the rectilinear vortex 
line. 

Another advantage of the quantum language in weak-turbulence problems, 
which we intensively employ in this paper, is that the collision term of the ki- 
netic equation immediately follows from the Golden Rule for the corresponding 
elementary processes. 



3 Pure Kelvin-wave cascade 

3.1 Qualitative analysis 

In the problem of low-temperature ST decay at the length scales smaller than 
the typical interline separation, the pure Kelvin-wave cascade on individual vor- 
tex lines is perhaps the most natural decay scenario one can think of. Indeed, 
the non-linear nature of Kelvin-wave dynamics should allow such a process in 
which the energy in the form of Kelvin waves is transferred from the long-wave 
energy-containing modes to the short wavelengths where the dissipation becomes 
efficient. Provided with a wide inertial range, one could also expect such a trans- 
fer of energy to be local in the wavenumber space, in which case the relaxation is 
likely to be due to a cascade. 

A subtlety that is immediately clear is the complete absence of kinetics in 
the leading approximation, which stems from the aforementioned integrability of 
the LIA, Eq. jilt . If it exists, the pure Kelvin-wave cascade must be entirely due 
to the non-local coupling between different vortex-line elements. Hence, being 
suppressed by the small parameter A " ' ^ 1 relative to the leading LIA dynam- 
ics, the pure KW cascade is an a priori weak phenomenon. This fact provides 
us with an important consistency check of the results: the leading contribution 
oc A — determined by the effective microscopic cutoff parameter a^, if one uses 
the pseudo-Hamiltonian ( 129b — must completely drop out of the effective kelvon 
scattering amplitude. 

The success of the Kolmogorov-type argumentation makes it very tempting to 
use a dimensional analysis of the dynamical equations of motion to immediately 
extract, e.g., energy spectra. Such a simple approach for Kelvin-wave turbulence, 
however, is doomed to failure. The reason is that here the scale invariance dictates 
the kinetics to be controlled by a purely geometrical diniensionless parameter = 
biJi, where bj^ is the typical KW amplitude at the wavenumber k. Thus, the solution 
could be only obtained from the corresponding kinetic equation, which reveals the 
proper combination of Uk and k the kinetics are due to. Deriving and solving the 
kinetic equation will be our main goal in the this section. 

Let us first determine the general structure of the collisional term. The key 
aspects that essentially fix this structure are (i) the conservation of the number 
of kelvons, discussed in Sec. [2] (ii) the dimensionality of the problem, and (iii) 
weakness of non-linearities. The first circumstance means that the kinetics are 
entirely due to kelvon elastic scattering since the elementary events of kelvon cre- 
ation/annihilation are prohibited. In other words, the only allowed kinetic channel 
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is energy-momentum exchange between kelvons. On the other hand, the energy- 
momentum conservation laws in one dimension make this exchange impossible 
in two-kelvon collisions: the process (^1,^2) ~> (^3 1^4) is only allowed if either 
(^3 —ki, ^4 = ^2); or (^3 = ^2i ^4 = ki) which does not lead to any kinetics. Thus 
any non-trivial kinetics can be only due to collisions of three or more kelvons. 

The third condition translates into the fact that the amplitudes of Kelvin waves 
in the pure cascade are necessarily small compared to their wavelengths, = 
bj^k ^ 1, at least at sufficiently large k. This is the central feature of this regime, 
which is due to the following. The cascade solution with the largest possible am- 
plitudes, a^. w 1 , corresponds to the regime driven by self-reconnections, as de- 
scribed in Sec.m the regime where the non-linear effects are negligible altogether. 
The weakness of the purely non-linear kinetics makes the amplitude at a scale k 
rise (due to the energy supplied from the larger length scales) until bi^ is of order 
l/k and the self-reconnection can happen causing the energy transfer to a smaller 
scale. If non-linear processes are appreciable, smaller amplitudes bj^ are sufficient 
to sustain the energy flux. Thus, the pure cascade spectrum, which due to the scale 
invariance has a general power-law form = bj^k oc kl^, must be constrained by 
j3 < 0. Apart from the marginal j3 = 0, the latter requirement guarantees that a the- 
ory built on tti <C 1 becomes asymptotically exact at hight wavenumbers. Corre- 
spondingly, self-reconnections, being exponentially dependant on ak, necessarily 
cease in the purely non-linear cascade. 

The smallness of the Kelvin-wave amplitudes in the particle language means 
that the many-kelvon collisions are rare events and we can confine ourselves to 
the leading three-kelvon processes in the kinetic term. 

These simple considerations already substantially limit the freedom in writing 
the collision term. In fact, the only missing ingredient, which needs to be cal- 
culated, is the wavenumber dependence of the effective kelvon scattering vertex. 
Postponing the final issue of finding this dependence until the next section, let 
us analyze the kinetic equation in a general form. Written in terms of averaged 
over the statistical ensemble kelvon occupation numbers ni^ = {a^j^aj^), the kinetic 
equation is given by 

^^=jf^i:^i^m-<ii)- (35) 

Here is the probability per unit time of the elementary three-kelvon scat- 

tering event (^1,^2,^3) {k4,ks,k(,), and the combinatorial factor compensates 
multiple counting of the same scattering event. The three-kelvon effective interac- 
tion Hamiltonian has the general form 

^^int= X! S{Ak)Vi2j alalalaj,a2ai, (36) 

k],...,k(, 

where the effective vertex V is symmetrized with respect to the corresponding 
momenta permutations, d{k) is understood as the discrete dj^ o, and Ak = ki + 
k2 + k^ — k4 — ks — k^. The probabilities Wj'*2 f thsn straightforwardly given by 
the Fermi Golden Rule applied to the Hamiltonian ( 136b : 

wt^:^=2n\i3l)'vtii\'jtip{Aoi)diAk), 
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/ui = nin2«3(n4 + l)(n5 + l)(«6 + l), (37) 
Am = 0)1 + 0)2 + 0)3-604-0)5-0)6. 

Here the combinatorial factor (3!)^ accounts for the addition of equivalent am- 
plitudes. We are of course interested in the classical-field limit of the Eq. ( 137b . 
which is obtained by taking n^. ^ 1 and retaining only the leading in terms. 
This procedure finally yields 

m=2l6n £ \vt:^:^\H{Aco)diAk) (/iii-/?!'' ) , (38) 

k2,...M(, 

/u,3 = «1«2«3(«4«5 +«4n6 + «5«6)- 

The kinetic equation ( l38l ) supports an energy cascade provided two con- 
ditions are met: (i) the kinetic time is getting progressively smaller (vanishes) in 
the limit of large wavenumbers, (ii) the collision term is local in the wavenumber 
space, that is the relevant scattering events are only those where all the kelvon 
momenta are of the same order of magnitude. In the following subsection, we 
make sure that both conditions are satisfied: the condition (i) can be checked by 
a dimensional estimate, provided (ii) is true. The condition (ii) is verified numeri- 
cally. Under these conditions one can establish the cascade spectrum by a simple 
dimensional analysis of the kinetic equation. 

The locality of the collision term implies that the integral in ( 1381) builds up 
around (^2, • ■ • ,k(,) ^ ki, which dramatically simplifies the kinetic equation: 

nk<^k^-\V\^-co^^^-k-^-nl, (39) 

where the factors go in the order of the appearance of corresponding terms in ([38] l. 
At ki ^ . . . ^ k(, ^ k we have |y| fc^ and 

fik oc o)-in|/t^+2^ . (40) 

The energy flux (per unit vortex line length), dk, at the momentum scale k is de- 
fined as 

Ok = L ' (41) 

k'<k 

implying the estimate 9^. ^ kfij^cOk- Combined with ( |4Q] i. this yields 0^ ^ «^^^+^^, 
and the cascade requirement that dj^ be actually fc-independent leads to the spec- 
trum 

n^oc (wt^^)=^^-(5+2v)/5 ^ (42) 

The value of the spectrum amplitude A in ( 1421) controls the energy flux that the 
cascade transports. The relation between 6 and A is 

d=CeK^pA^. (43) 

The dimensionless coefficient in this formula can be obtained, e. g., b y a straight- 
forward numerical calculation, as it was done by the authors in Ref. 11311 . However, 
the calculation of Ref. fl3ll in view of a large relative error allowed to obtain only 
the order of magnitude of Cg, which, as we explain in the next subsection, can 
be further questioned due to an erroneous omission of an order-one term in the 
effective vertex. 
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3.2 Quantitative analysis 

The central problem of this section is the derivation of the effective vertex V^2 3^ 

which defines the effective kelvon interaction Hamiltonian ( 136b {V is obtained 
from V by symmetrization with respect to corresponding momenta permutations). 
We start with the pseudo-Hamiltonian ( 129b . Our fundamental requirement that the 
amplitude of KW turbulence is small as compared to the wavelength, ^ 1 , is 
formulated by Eq. ( l30l ). This allows us to expand ( |29l ) in powers of oc{zi,Z2) ^ 
1 : H = Eo + Hq + Hi + H2 + ■ ■ ■ {Eq is just a number and will be ignored). The 
term Hq is given by Eq. i3li . As we demonstrated in the previous subsection, it 
describes the linear properties of the Kelvin waves, in particular it determines the 
Kelvin-wave dispersion law, Eq. ( 1331) . The higher-order terms are responsible for 
interactions between kelvons. The terms that will prove relevant are 

Hi^:^ ( ^^[3a4-4a2Rew'*(zi)w'(z2)], (44) 
3271 J \z\~Z2\ 

|zi-?2|>a« 

and 

^2 = 7^ / ^^i^[6a4Rew'*(zi)w'(z2)-5a'']. (45) 

0471 J \Z\—Z2\ 

\z\-Zl\>a* 

As already demonstrated, the leading elementary process in our case is the 
three-kelvon scattering and the processes involving four and more kelvons are 
much weaker due to the inequality ( l30l ). The effective vertex, Vj 2 3 , for the three- 
kelvon scattering [subscripts (superscripts) stand for the initial (final) momenta] 
consists of two different parts. The first part is due to the terms generated by the 
two-kelvon vertex, A (corresponding to the Hamiltonian /fi ) in the second order 
of perturbation theory. [All these terms are similar to each other; we explicitly 
specify just one of them: A^2G{a>i ,k-i)A-i'T^. Here G{(0,k) = l/(co — (Ok) is the 
free-kelvon propagator, CO7 = 0)1 + ft>2 ~ ^4, h = ki + k2 — ^4.] The second part 
of the vertex V is the bare three-kelvon vertex, B, associated with the Hamiltonian 
7/2. The expression for the effective vertex is depicted diagrammatically in Fig.|2l 

The explicit expressions for the bare vertices directly follow from (144b and 
( |45] | after the Fourier transform w{z) = L'^I^Y.k'^k^''^'"- 

A = {6D-E)/%n, (46) 
Elj = j\dx/x'){kM['] + [,]-['']-i]) + 

^3fci(rV[j-n-[2])+wrV[2]-n-ffl) 
+^4M[V[2]-n-[2])}^ 
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Fig. 2 The effective vertex of the three-kelvon scattering. The vertices A and B come from the 
terms Hi i44\ and H2 <45t and are given by <46t and <47t respectively. G{co,k) is the free kelvon 
propagator. 



5=(3P-5e)/47r, (47) 

eti3=/VA^){i-r]-[j+B-[vr] 



Here [ ] 's denote similar looking cosine functions, [J =cos^ix, [J =cos(fc4 — 

[j^] = cos(^4 + ^5 — ki)x, [jj] = cos(fc4 + ^5 — fci — fe)^^, and so forth. 
Since the microscopic physics at the scales ^ uq is quite complicated and usu- 
ally not known precisely in strongly correlated systems like ^He, it is important to 
asses the systematic error of the theory in the case when a* can not be determined 
accurately. For this purpose, as well as for the convenience of the numerical anal- 
ysis described below, we process the integrals ( 146b , ( 147b as follows. We introduce 
a characteristic wavelength X = X{k\,k2,kT,,kn,k^,k(,), the particular functional 
form being not important. Then, from each of the integrals ( 146b , ( 147b , we sub- 
tract the corresponding LIA contributions, which are easily obtained by the same 
procedure that led to Eqs. ( 1461 ), ( 147b from the local Hamiltonian ( 1241 ) defined by 



j3;L = ('c/27r)A;L, A^=\n{l/a,). 



(48) 
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As a result of the subtraction, up to system-specific terms ^ to be ne- 

glected below, we get a* -independent convergent integrals. Thereby we arrive at a 
convenient decomposition: A = A^A^^^ +A^^\ B = A^B^'^^ +B^^\ where A'*') and 
5(0) 

are known analytically from LIA, while A'^^' and B^^^ are easily calculated 
numerically (see below for some useful details). On the technical side, the de- 
composition solves the problem of handling the logarithmic divergency, whereas 
physically, it clarifies the form of the dependence on the vortex-core details, which 
turn out to be entirely enclosed by LIA. 

To proceed with estimating possible systematic errors, we introduce the ex- 
pansion of the propagator G in the powers of inverse Ax : 

G = A;^1[g(") +A;[1g(>) + e{A^% (49) 

For the vertex V = AGA + B (see Fig.|2]l, we thus have 

V=AxV'^^'> +y(') + V(2)/A;i + &{A^'^), (50) 

with 

y{0) ^a^^^g'^U'") 
yd) = 2aWgWa(') +aWg(')a(«) 

yd) = aWg(2)a(«) +2aWg('U(') +a(')g(''U(i). 

Note that by construction all the quantities A^'^ fi^'', g''\ and V are a^, -independent, 
since the dependence on — up to the neglected terms ^ (a^./X)'^ — comes exclu- 
sively through A^. The term A;^y'''' is precisely the effective vertex that follows 
from the LIA Hamiltonian (|24] |. and thus it necessarily obeys V = 0. The lead- 
ing contribution to V is the a* -independent V^^\ Thus, the short-range physics 
enters the answer only as a small in 1 /A^ correction. 

Now we are in a position to specify the systematic error due to uncertainties 
in microscopic details. In the case when it is possible to calibrate the pseudo- 
Hamiltonian, that is to find (analytically or experimentally, see Sec.|2]l the accurate 
value of fl*, the systematic error is of order Otherwise, we are forced to 

set fl* ^ ao, meaning that the parameter A;^ ^ 1 is only known up to 8Ai ^ 1, 
which for the systematic uncertainty in V gives 

y(2)5A;LMi ~ IMa- thus 
arrive at an important conclusion that even if the microscopic details, such as the 
vortex-core shape, are completely unknown, by naively setting A^ = ln(/l/flo) we 
are paying by a relative error of only 1 /A|. 

A comment is in order here on the technical issue of a convenient handling 
of the integrals, which we do numerically, upon subtracting the logarithmic sin- 
gularity. The whole procedure is almost identical to the one that led us to the 
kelvon dispersion, Eqs. (I33b -( l34b . The new aspect is the dependence on four/six 
momenta, rendering a direct tabulation of integrals computationally expensive. 
The trick is to split a four/six-parametric integral into a sum of single-parametric 
integrals. A minor technical problem comes from the power-law divergence at 
X — > of each separate single-parametric integral in Eqs. ( 146b , ( 147b . The problem 
is readily solved by introducing power-law counter-terms — the net contribution 
of which is identically equal to zero — rendering each individual single -parametric 



18 



integral convergent. As an illustration, we present the final result for the regular- 
ized (by subtracting the LIA terms and introducing the counter-terms) integral D 
of Eq. ( 146b , which we denote with D. With the re-scaled (by x — > Xx) integration 
variable, the expression reads 

= x-'j\dxix% 1 - [J - u - - + a + r^] + a) 
j\dx/x'){\^,]] + n + [fl] - u - [[,]] - [f]] - n). (51) 

Here [ ] 's denote the same cosines as previously, but with an extra factor X in 
the argument due to re-scaled x. The symbol [[ ]] means that for the cosine, we 
subtract first few — first three in the case of D — terms of its Taylor expansion to 
render corresponding single-parametric integral convergent (and thus individually 
tabulatable). 

From Eqs. (146b , ( 147b it is straightforward to obtain the scaling of the effective 
vertex V with the momenta — at k\ ^ . . . ^ k(,^ k^Ne. have \ V\ ^ k^. Thus, v = 6 in 
Eq. (142b and the pure Kelvin-wave cascade spectrum (restoring all the dimensional 
coefficients) is 

{wl^,) = ^=Ak-"/\ (52) 

This result was corroborated in a direct numerical simulation by the authors 
fl4ll . where the spectrum ( 152b was resolved with a high accuracy allowing us to 
distinguish it from n oc k~^ suggested in an earlier simulation by Vinen et al. 
|[l2ll . The high precision required to distinguish the close exponents, both in terms 
of the cascade inertial range extent and low noise, was achieved using a special 
numerical scheme developed to reduce the complexity of the non-local model ( 1291 ) 
to that of an effectively local one at an expense of a controllable systematic error. 
We refer to Ref. fl3ll for more details. 

Let us also express the spectrum in terms of the typical geometrical amplitude, 
b]^, of the KW turbulence at the wavevector ^ k. By the definition of the field w{z) 
we have; ^ L~^Y,q-^k{wlwk) k (wlwk). Hence, using Eq. ( |43] |. we obtain 

bk-{0/K^p)^f^h-^f\ (53) 

One can convert ( |52] i into the curvature spectrum. For the curvature c{Q = d^s/dl^^ 
[where s(i^) is the radius-vector of the curve as a function of the arc length the 
spectrum is defined as Fourier decomposition of the integral Ic — J \c{Q\-^d^. The 
smallness of a, Eq. dSOl l, allows one to write Ic ^ J dziw""^ {z)w" (z)) = Y^k^^'^'^k 
^i^fc^/^, arriving thus at the exponent 3/5. 

So far we were heavily relying on the assurnption of locality of the kinetic pro- 
cesses in the wavenumber space. We checked |[l3n the validity of this assumption 
by a numerical analysis of the kinetic equation thereby (i) making sure that the 
collision term of the kinetic equation is local and (ii) estimating the value of the 
dimensionless coefficient in (143b (see, however, below). The analysis is based on 
the following idea ll42h . Consider a power-law distribution of occupation numbers, 
ri]^ = Ajk^ , with the exponent j3 arbitrarily close, but not equal, to the cascade ex- 
ponent j3o = 17/5. Substitute this distribution in the collision term of the kinetic 
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equation — right-hand side of ( 138b . Given the scale invariance of the power-law 
distribution, the following alternative takes place. Case (1): collision integral con- 
verges for j3's close to /3o, and, in accordance with a straightforward dimensional 
analysis, is equal to 

Con{[nk=A/kP],k)=C{l5)A^0),;^/k^P-^^ . (54) 

Here C(j3) is a dimensionless function of /3, such that C(/3o) = since the cascade 
is a steady-state solution. Case (2): collision integral diverges for j3 close to j3o. 
The case (2) means that the collision term is non-local and the whole analysis 
in terms of the Kolmogorov-like cascade is irrelevant. Fortunately, our numerics 
show that we are dealing with the case (1). Substituting ( 1541 ) for % in ( 141b . we 
obtain the expression 0<.(j3) = {Ay2n)k"-^P C(j3)/(17 - 5j3). Taking the limit 
P j3o, we arrive at the fc-independent flux 6 = -C' {l5o)A^ /IQn. 

In Ref. |(13i1 , we used this formula to obtain the coefficient in ( 1431 ) by calculat- 
ing C(j3) and finding its derivative C'(j3o) = — IOttCq. We simulated the collision 
integral by Monte Carlo method. The integrals ( 146b , ( 147b . were calculated numer- 
ically. However, a mistake was made at the level of combining these integrals into 
the effective vertex V. Due to the cancelation of the leading logarithmic terms 
one has to keep an accurate account of the sub-logarithmic corrections, includ- 
ing those in the kelvon propagator G{(0,k) = l/(ft) — 0)^) as well, as we already 
demonstrated above. Unfortunately, we failed to appreciate a simple fact that the 
constant Co appearing in the kelvon dispersion (133b results in an order-one con- 
tribution to the effective vertex V, and neglected this constant in the propagator. 
Since the coefficient Ce ~ 10^^ was obtained in Ref. fl3ll only as an order of 
magnitude estimate (with an error ^15% due to a slowing down of the high-order 
numerical integration) the mistake might not significantly effect this result, but 
still makes it questionable. A more accurate determination of Cg is necessary. 



4 Self-reconnection driven cascade 

4.1 Absence of Feynman's cascade 

For quite a long period of time it has been generally accepted that at T = the 
scenario of decay of superfluid turbulence is the one proposed by Feynman 13^ . 
Namely, that vortex lines first decay into vortex rings, each of the rings then inde- 
pendently self -reconnects, producing smaller rings, each of the smaller rings self- 
reconnects to produce even smaller rings, and so forth. In this respect it is very 
characteristic that even after the absence of pure Feynman's cascade was directly 
shown, and the scenario of Kelvin-wave cascade driven by local self-crossings was 
proposed isj], simulations of the low-temperature decay within LIA performed in 
Ref. i] were still interpreted in terms of Feynman's cascade. 

The problem with Feynman's cascade is that it is inconsistent with simultane- 
ous conservation of energy and momentum @l- This is clear from the mere fact 
that energy scales as the length of the ring (up to logarithmic corrections) while the 
momentum scales as the length squared, so that a decay into arbitrarily small rings 
with the net line length conserved would result in vanishing total momentum. 
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4.2 Generation of Kelvin waves in the process of line reconnection 

With a slight modification of the time dependence of the phase, the standard 
self-similar solution of the linear Schrodinger's equation — its Green's function — 
applies also to the non-linear equation ( 1161) . Indeed, the absolute value of the solu- 
tion is spatially homogeneous, so that the non-linearity is equivalent to a spatially 
homogeneous time dependent external potential that can be immediately absorbed 
into the phase as an extra time-dependent term. The result is 



(11 a2 



(55) 



In accordance with ( 1151 ), this solution implies 

^=A/Vt, (56) 
T = ^/2f. (57) 

The physical meaning of the solution ( |56I )-(I57||, first revealed by Buttke ll43h . is 
the relaxation of the vortex angle, the value of which is controlled by the param- 
eter A. This solution gives an accurate description (within LI A) of relaxation of 
two vortex lines after their reconnection, as long as the curvature in the relax- 
ation region remains much larger than the curvatures of the two lines away from 
the region, so that the distant parts of the two lines can be approximately treated 
as straight lines. By dimensional argument, the self-similarity regime should be 
achieved very rapidly upon the reconnection, with the velocity of propagation of 
the fastest Kelvin waves (with wave vectors flo)- 

With the solution (I56ll-(I57|) one can explicitly see that reconnections lift the 
integrability constraints. At ^ oo we have (p„ ^ meaning that a single recon- 
nection renders all the integrals /„ dlSb divergent. 

Let us look at the asymptotic form of the solution (I56l)-(l57]) in Cartesian co- 
ordinates, taking the direction of the z-axis along one of the two lines, z +°° 
corresponding to the asymptotic limit: 

x{z,t)+iy{z,t) = {4Aty^/z^y''/^', z->Vt- (58) 

Equation (158b reveals a helical Kelvin-wave structure moving away from the re- 
connection region. It is important that while arbitrarily small wavelengths are 
present in the solution ( 158b , the integral for the total line lengths comes from 
the largest length scale z ^ \ft- This provides a support for two crucial points of 
the scenario of self-reconnectionsdriven Kelvin-wave cascade: (i) Reconnections 
push Kelvin waves to smaller lengthscales, (ii) apart from higher-order correc- 
tions, the line length associated with the curvature radius R cannot go directly to 
the scales of curvature radius much smaller than R (locality of the cascade in the 
wavenumber space). 

When the angle between the two lines is close to tt, the reconnection leads to 
a production of vortex rings ll43ll . This process does not introduce a new cascade 
channel. Qualitatively it is very similar to the pure production of Kelvin waves 
( 158b . because, up to higher-order corrections, the line length is being transferred 
to the rings of the radii of the order of the curvature radii of reconnecting lines, 
i.e. to adjacent scales in the kelvon wavenumber space. 



21 



4.3 Fractalization of lines. Kelvin-wave spectrum 

In this subsection we render the analysis of Ref. lHj] of the fractalization of a vortex 
line by the self-crossings driven Kelvin-wave cascade. We start with introducing 
the crucial notion of a smoothed line length, ^(Ai), which is the length of a 
(fractalized) line upon smoothing out all the structures of the length scales smaller 
than Xi . Corresponding mathematical expression reads 

In^(Ai) - In^(Ao) + {b^/XfdX /k, (59) 

where is the characteristic amplitude of the Kelvin-wave structure at the wave- 
length A, while Ao is, generally speaking, any fixed wavelength scale significantly 
larger than X\ . In particular, if Ao is the largest wavelength of the problem, then 
^(Aq) is the length of maximally smoothed vortex line. With J^{X) we can esti- 
mate the number of local self-crossings at the scale A per unit time as 

N^^^^^1{X)C0^, (60) 

where the factor ^(A)/A gives the number of statistically independent pieces of 
the line (at the wavelength scale A typical correlation length is ^ A), the factor 
i2(A) is the probability to have a large enough amplitude ^ A to produce a self- 
crossing within a given element of the line with the length ^ A, and (Ox is the 
kelvon frequency playing the role of inverse correlation time. 

To estimate X2(A), we can rely on the theoretical limit of bx ^ A, in which 
kelvons are independent harmonic modes and, correspondingly, the statistics of 
fluctuations of the amplitude is Gaussian. This readily yields 

^2(A)-(^;l/A)e-(^/''^)^ (61) 

We do not introduce any dimensionless factor in the Gaussian exponent in view of 
the freedom of defining bx up to a factor of order unity. 

For any cascade the fundamental notion is the flux of corresponding conserved 
quantity. In our case, it is the flux of the vortex line length. The fractalization of the 
lines introduces certain subtleties. In contrast to a standard cascade in which the 
integral for the conserved quantity comes from a single length scale, in our case 
the line length is spread over all scales of distance of the inertial range. Moreover, 
different wavelength scales are not entirely independent in the sense that with the 
fractalized lines the short-range structures with their energy are slaved by the long- 
wave structures. So that the line length coming from long waves to shorter ones 
is essentially carried by short-ranged structures slaved by the long-wave modes. 
The crucial observation now is that with respect to larger wavelengths the slaved 
shorter wavelengths play a rather passive role in energy balance, since short-wave 
contribution to line length carried by the longer waves is just proportional to the 
smoothed line length. Correspondingly, one can speak of the smoothed line length 
flux, Q(Ai ), where Ai is the smoothing parameter of Eq. ( |59t . By definition of the 
cascade, the quantity Q is one and the same for any wavelength scale A, as long 
as A ^ Ai . 
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Having fixed some wavelength scale A, we note that the self-reconnections 
scenario implies 

e(Ai)^7V;i^;i(Ai), (62) 

where S^^ {X\ ) is the length of a /li -smoothed circle of the radius ^ A. In a direct 
analogy with ( 1591) . we have 

ln^;i(Ai) -lnA+ / {b}_, /X'fdX' /X' . (63) 

From ( 1631 ) and ( |59l l the follows a useful relation 

J^{X)^.x{h)-^-^{h). (64) 

WithA^;L (llQll , ^2(A) (|6T|, the estimate CO;^ i^/'^^' and the relation (|64l), equation 
( |62] | yields 

e(Ai)/^(Ai) ^ (/3^;i/A3)e-(V''A)'. (65) 

The left-hand side of this relation is a function of X\ , while the right-hand side is 
a function of A, meaning that both sides are actually constants. For the right-hand 
side this implies 

(Since up to logarithmic corrections we have ^ A, we do not distinguish be- 
tween \nbx and InA.) Expression ( 166b yields the Kelvin-wave cascade spectrum 
in terms of the characteristic amplitude bx. As far as the function ^(A) is con- 
cerned, from ( 166b and ( [59l ) we find 



if(A)=if(Afl) 



^Vln^ 
Ao 



V 

(67) 



where V is a constant of order unity the particular value of which cannot be estab- 
lished by our order-of-magnitude analysis. 

It is clear from Eq. ( |53l l that no matter how large the energy flux (per unit 
vortex-line length) 9 °^ Q transported by the self-crossings regime is, at suffi- 
ciently high wavenumbers the pure Kelvin-wave cascade will be capable of sup- 
porting it. As soon as the purely non-linear kinetics become appreciable, the am- 
plitudes of Kelvin waves must become smaller, which in view of Eq. ( 1611) , inhibits 
the reconnections. Thus, the self-crossings-driven regime will inevitably be re- 
placed by the pure Kelvin-wave cascade at some scale A*. Note, however, that 
due to a small difference between the Kelvin-wave spectra in the two regimes, 
Eqs. ( |53l l, ( |66l l, the crossover between them is likely to be extended in the wavenum- 
ber space. A rough estimate of the scale A* can be obtained by setting b^^k"^ ^ 
A* in Eq. ( 153b . The result clearly depends on the cascade energy flux 9, which 
is specific to the physics at the energy-containing scale. In quasi-classical turbu- 
lence, 9 is related to the Kolmogorov energy flux, in which case an estimate for A* 
will be obtained in Sec.|5] In the case of a non-structured tangle, the energy flux 
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is formed by reconnections at the scale of interline separation /q, ©ns ^^P^^/^o^ 
with A = ln(/o/a*), which gives 

^* ^ h/-^ (non — structured tangles). (68) 

According to this rough estimate, the inertial range for the regime driven by local 
self-crossings is only about A/ 10 decades, which for realistic values of A could 
turn out to be an insignificant range without a distinct spectral signature. Thus, it 
would be crucial to quantify the crossover between the regimes by a direct numeric 
simulation. 



4.4 Qualitative Hamiltonian model 

The Hamiltonian description ( 1191 ) implies single-valuedness of the function w(z). 
Amazingly, the question of what happens to the mathematical solution of Eq. ( 1191 ) 
when the physical function w{z) becomes non-single valued turns out to be very 
relevant to the theory of self-crossings driven cascade. Clearly enough, the math- 
ematical solution has to develop a certain singularity, and, strictly speaking, be- 
come ill defined afterwards. However, if one introduces a discretized analog of the 
problem described by the Hamiltonian fsj] (we confine ourselves to LIA), 



iwn = 3—, H = y + \w„+i - w„|2, (69) 

the Hamiltonian dynamics remains well-defined at any time moment. Numeri- 
cal analysis of this model @] shown that that the above-mentioned singularity 
evolves into a finite jump between the values of w j and Wy+i, at a certain 7. The 
jump exists for a certain time, and then relaxes, the relaxation process being qual- 
itatively similar to the process of vortex angle evolution in the sense that Kelvin 
waves are being emitted, while the integrability constraint is naturally lifted by 
non-smothness of the function. With the precise geometric meaning of the Hamil- 
tonian (|69l l — the length of the broken line defined by the points {w„} — we realize 
that the dynamics governed by it should be qualitatively equivalent to the self- 
reconnection induced cascade, leading to the fractalization of the line necessary 
to support jumps of arbitrarily small amplitude, the amplitude of the jump playing 
qualitatively the same role as the radius of the ring in the self-crossings driven 
cascade. And that is precisely what has been revealed by numeric simulation of 
the model W)\ in Ref. IH]. The simulations also revealed the spectrum ^ k, 
consistent with Eq. (166b , the logarithmic factor going beyond numeric resolution. 



5 Crossover from Richardson-Kolmogorov to Kelvin-wave cascade 

5.1 Quasi-classical tangles at T=0: crossover to the quantized regime 

Even at absolute zero temperature the superfluid dynamics supports a turbulent 
regime, which under certain conditions is indistinguishable from classical turbu- 
lence I4i]. That may seem quite surprising since the only degrees of freedom in 
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a superfluid at T = are quantized vortex lines, which are singular topologi- 
cal objects and thus are very different from the classical eddies responsible for 
turbulence in classical ideal incompressible fluids. Nonetheless, quantized vortex 
lines possess a mechanism that allows them to mimic classical vorticity — it is well 
known jj]] that macroscopic velocity profile of a rapidly rotated superfluid mim- 
ics solid-body rotation, which is accomplished by formation of a dense array of 
vortex lines aligned along the rotation axis. The basis of this mechanism is the 
strong coupling between the vortex lines in the dense array, which makes such 
a bundle behave as a single coherent classical object. Therefore, by essentially 
classical turbulence generation methods (i.e. "stirring") one can produce vorticity 
in the course-grained up to length scales larger than the typical interline separa- 
tion If) superfluid velocity field, indistinguishable from that of a normal fluid, the 
underlying vortex tangle being organized in polarized "bundles" of vortex lines. 

Over the last decade, experimental observations of the classical behavior ex- 
hibited by superfluids [29, 30, 31, 18, 16, 19, 24] have largely led to a renaissance 
of general interest in superfluid turbulence. Perhaps the most attractive feature of 
this quasi-classical turbulence is that, unlike counterflow turbulence, it in principle 
allows generation and probing at temperatures close to absolute zero. [It is only 
very recently that a unique technique of non-structured tangle (a la counterflow 
turbulence) production at very low temperatures was developed by the Manch- 
ester group iU].] With recent technological advances this opens an intriguing 
possibility of studying such essentially low-temperature phenomena as, e.g., the 
Kelvin-wave cascades. 

However, it was recently realized ll20ll that, at T = 0, the question of how 
the quasi-classical vortex tangle looks like when one zooms in down to scales of 
order Iq, where the vorticity is essentially discrete, is quite a puzzling problem. 
The only fact that is immediately clear is that near the scales Aph <C Iq where 
the dissipation due to the sound radiation takes place, the energy flux must be 
transported by the pure Kelvin-wave cascade. What happens in the intermediate 
regime between the classical Kolmogorov cascade of eddies and the Kelvin-wave 
cascade on individual vortex lines is the subject of this section. 

In their scenario, L'vov, Nazarenko, and Rudenko |2QJ, noted that a simple pic- 
ture in which the pure Kelvin-wave cascade supersedes the Kolmogorov regime 
at the scale Iq is not possible. The difficulty is due to the fact that at this scale the 
pure Kelvin-wave cascade is unable to sustain the Kolmogorov energy flux (per 
unit mass of the fluid) £. Correspondingly, L'vov, Nazarenko, and Rudenko put 
forward an idea of bottleneck accumulation of energy at the classical scales adja- 
cent to Iq. The accumulation of energy in the form of a thermalized distribution of 
quasi-classical vorticity was suggested to be necessary to raise the level of turbu- 
lence to a value at which the pure Kelvin-wave cascade becomes efficient. Note, 
however, that the concept of bottleneck at a given scale fundamentally relies on 
the absence of any efficient energy transport mechanism at this scale. Due to this 
fact, under the conditions of Ref. |[20ll . vortex-line reconnections play a fatal role 
for the bottleneck scenario. Indeed, if the Kolmogorov cascade can reach the scale 
Iq as assumed in Ref. ll20ll . the typical vortex-line curvature at this scale is of the 
order Iq meaning that the vortex-line reconnections must happen due to the tangle 
geometry, which makes them an alternative energy transport channel to the pure 
Kelvin-wave cascade. Estimating the energy flux £rec processed by the reconnec- 
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tions at this scale we see that Srec/s ^ 1, i.e. we are actually dealing with 

an "anti-bottleneck" — the reconnections transport an energy flux much larger than 
the one supplied from the larger scales. The anti-bottleneck is of course forbidden 
by the energy conservation, so we are forced to conclude that the transformation 
of the classical regime must happen already before the scale Iq is reached, and the 
reconnections play a crucial role in this process. 

Let us first analyze the tangle of vortex lines at T = at the length scales 
much larger than the typical interline separation under the condition of a developed 
Kolmogorov cascade. Although understood intuitively, the ability of the system of 
quantized vortex lines to mimic the classical Kolmogorov cascade is not trivial. 
We start by a rigorous demonstration of this fact. 

At r = 0, vortex lines are the only degrees of freedom and their dynam- 
ics are completely captured by the Biot-Savart equation The latter can be 
rewritten in classical terms of vorticity w = curlv in the momentum space, Wk = 
/ w(r) exp[— /k • r] d^r = kJ exp[— /k • s] ds. The result is identical to the vorticity 
equation for a normal ideal incompressible fluid: 



In view of Eq. ( 170b we can formulate the conditions, under which the vortex tan- 
gle must be automatically equivalent to classical ideal-incompressible-fluid turbu- 
lence: (i) the energy must be concentrated at a sufficiently small wavenumber scale 
^en ^0 ^ ^nd (ii) the decay scenario must be local in the momentum space, so 
that the quantized nature of vorticity is irrelevant for the long-wavelength behav- 
ior. These conditions are not restrictive: (i) is automatically satisfied if turbulence 
is generated by classical means due to the large compared to k values of the veloc- 
ity circulation, and (ii) is necessary for the existence of the Kolmogorov cascade 
in classical fluids as well. 

Note that the circulation quantum K completely drops out of the vorticity equa- 
tion. This is a manifestation of the known fact that the superfluid hydrodynamics 
is completely described by the classical Euler equation with respect to which the 
quantization of circulation is nothing but an imposed initial condition enforced by 
quantum mechanics — due to the Kelvin theorem, once preformed the velocity cir- 
culation is a constant of motion. Since the dynamics of each individual vortex line 
are controlled by the circulation quantum k, the independence of Eq. ( |70b on K 
leads to an important conclusion that any large-scale (classical) motion necessarily 
implies strong coupling of the underlying vortex lines and that the crossover to the 
quantized regime is due to the self-induced motion of the vortex lines starting to 
dominate over the inter-line coupling. To obtain the corresponding crossover scale 
ro, let us formally decompose the integral (|4]l into the self-induced part, v^'(s), for 
which the integration is restricted to the vortex line containing the element s, and 
the remaining contribution induced by all the other lines v^(s). 



By the definition of ro, at length scales the turbulence mimics classical 

vorticity taking on the form of a dense coherently moving array of vortex lines 




(70) 



v(s)=vS'(s) + v'(s). 



(71) 
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bent at a curvature radius of order r. The velocity field of this configuration obeys 
the Kolmogorov law 

Vr^{er)'l\ ry^ro, (72) 
where £ is the energy flux per unit mass of the fluid formed at the energy-containing 
eddies and transferred by the cascade. Here and below the subscript r means typ- 
ical variation of a field over a distance ^ r. On the other hand, the value of Vr 
is fixed by the quantization of velocity circulation around a contour of radius r, 
namely Vyr ^ Kriyr^, where n,- is the areal density of vortex lines responsible for 
the vorticity at the scale r. Note that scale invariance requires that on top of vor- 
ticity at the scale r there be a fine structure of vortex bundles of smaller sizes, so 
that, mathematically, n, r^ is the difference between large numbers of vortex lines 
crossing the area of the contour r in opposite directions. The quantity is related 
to the flux by 

1/3 

, r > ro . (73) 



The underlying dynamics of a single vortex line in the bundle is governed by 
v]. and vf^. While by its definition ^ v,-, which is given by Eq. ( 172b . the self- 
induced part is determined by the curvature radius r of the vortex line according 
to the LIA, Eq. (dB, 



vf-A,-, (74) 
r 

where = ln(r/a*). Here and throughout this section, the logarithmic accu- 
racy will be sufficient for the analysis, and, correspondingly, we replace Ar with 
A = ln(Zo /oq). At length scales where vl ^ vf^, the vortex lines in the bundle move 
coherently with the same velocity ^ v].. However, at the scale rg (A^)C^/e)'/^, 
the self-induced motion of the vortex line becomes comparable to the collective 
motion, vf^ ^ v].. At this scale, individual vortex lines start to behave indepen- 
dently from each other and thus ro gives the lower cutoff of the inertial region of 
the Kolmogorov spectrum (172b . 

Since ro is the size of the smallest classical eddies, the areal density of the 
vortex lines at this scale is given by the typical interline separation, n,Q ^ l/^o- 
other words, vortex bundles at the scale ro consist of almost parallel vortex lines 
separated by the spacing Zq. With Eq. (173b , we arrive at 

(75) 

/o-(AjcVe)'/^ (76) 
The crossover can be also understood in slightly more visual terms. Let us 
introduce an effective number of vortex lines A^,- in a bundle of size r. This number 
is obtained as an algebraic sum of the number of lines going through the bundle 
cross-section in opposite directions and is given K. In 

view of Eqs. (|72] |. ( |74] |. the number of lines in a bundle relative to A determines 
whether it behaves as a classical eddy or a set of independent vortex lines: A^,- ^ 
A means that the coupling between the vortex lines dominates resulting in the 
crossover when A^, „ ^ A . Thus, in the theoretical limit of A ^ 1 the bundles still 
contain a large number of vortex lines at the scale where the classical regime 
breaks down. 

In the next subsection we shall introduce the cascade mechanism that super- 
sedes the Kolmogorov cascade of eddies at the scales immediately adjacent to ro. 



K 
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5.2 Reconnections of bundles 

The analysis presented in this and the following section relies on the fact that the 
number of vortex lines in a bundle at the crossover scale is large, A^^q ^ 1 , which, 
in view of Eq. ( 175b , is guaranteed in the limit of large A . 

At the scale ro, turbulence consists of randomly oriented vortex-line bundles 
of size ro formed by the classical regime. The length ro plays the role of a cor- 
relation radius in the sense that relative orientation of two vortex lines (with the 
short-wavelength structure smoothed out) becomes uncorrelated only if they are 
a distance > ro apart. On the other hand, the crossover to the quantized regime 
means that each line starts moving according to its geometric shape, as prescribed 
by Eq. dl lb . Therefore, reconnections, at least between separate bundles, are in- 
evitable and, as we show below, capable of sustaining the flux e. 

The quantity that will play an important role in the analysis is the energy trans- 
ferred to a lower scale after one reconnection of vortex lines at the scale k~^, 
which, following Ref. js}], can be written as 

e,^f{r)ApK'k-'. (77) 

Here, /(y) is a dimensionless function of the angle y at which the vortex lines 
cross (y = corresponds to parallel lines). Its asymptotic form is 

f{y)-y', y«i- (78) 

Although at the scale ro there is already no coupling between vortex lines to 
stabilize the bundles, they should still move coherently on the time scale of their 
turnover time since, by the definition of the bundle size rg, geometry of neigh- 
boring lines at this scale is essentially the same over distances < ro. On the other 
hand, during about one turnover the bundle must cross a neighboring bundle and 
reconnect providing a mechanism of energy transfer to the lower scales. It is pos- 
sible, however, that vortex lines within the bundle reconnect. One can show that 
such processes can not lead to any significant redistribution of energy at the scale 
ro (but they will play an important role at smaller scales) and thus to a deformation 
of the bundle at this scale because they happen at small angles so that the energy 
( 177b is too small. Indeed, the dimensional upper bound on the rate at which two 
lines at a distance / <C ro can cross each other is, from Eq. (II lb . A K/r^l, while 
the actual value should be much smaller due to the strong correlations between 
line geometries. Taking into account that the number of lines in the bundle is 
{fo/lo)^ and that l/ro, the contribution to the energy flux from these processes 
is bounded by (//ro)e. Only when I ^ ro the reconnections become important, 
which are the reconnections between the whole bundles of size ro. 

Crossing of the bundles results in reconnections between all their vortex lines 
and Kelvin waves with a smaller but adjacent wavelength A are generated. This 
picture of bundle crossing was recently corroborated by direct numerical simula- 
tions 125] . The coherence of the initial bundles implies that the waves on different 
vortex lines of the same bundle must be generated coherently. Thus, at the scale 
k^^ < ro, adjacent vortex lines should still be almost parallel with Kelvin waves 
on them of the wavenumber k and amplitudes Z?^ — vortex lines at the scale k~^ 
also form bundles. Similarly, these bundles can reconnect transporting the energy 
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to a lower scale, where similar bundle reconnections happen, and so on down the 
scales until the bundle size is comparable to the interline distance and the self- 
similar regime is cut off. 

An important ingredient of this scenario characterizing the bundles at a scale 

k is the correlation radius rl_ of the vortex-line geometry (with the short- 
wavelength structure smoothed out up to fc" ' ) in the transverse to the bundle di- 
rection. The value of determines the distance over which neighboring vortex 
lines can be considered as parallel, i.e. it gives the transverse size of the bundle. 
This size is due to a finite time required for the Kelvin-wave amplitude at the scale 
^ to build up (after a reconnection at a larger scale ' > ^2 ^ ), which is of order 
of the wave turnover time, Tj^^^^ \/ kA k\ . The Kelvin waves that were generated 
within the time ^ Tj^ are coherent. In other words, as a reconnection of two lines 
happens at time / = T^j at the scale k^^ > ^2 V the lines that reconnected in the 
same bundle at f = have already-developed waves, which can not be coherent 
with the waves about to be generated at f > . Thus, the distance traveled by the 
bundle at the scale k^ ' over the time ^ T^^ determines the orientational correla- 
tion radius at the scale k2, = r]^^ {k\) = bkik\/k\. Since the scales k\ and ki 
are actually adjacent (i.e. different only by a factor of order unity), we finally get 

The spectrum of Kelvin waves b]^ in this regime can be obtained from the 
condition Ek = £, where Ek is the energy flux per unit mass transported by the 
reconnections at the scale k~^, 

h-{klp[rt^]^)Nke^T^\ (79) 

Here, we take into account that the correlation volume of the reconnection is 
[r'^\'^ /k, ^ {bk/lo)^ is the number of vortex lines participating in the recon- 
nection, and T^~' KAk^ is the rate at which the bundles cross. Physically, b^ 
determines the typical crossing angle, bf^k, thereby controlling the energy lost 
in one reconnection. From Eq. ( 1791 ), the spectrum of Kelvin waves in the bundle- 
crossing regime has the form 

bk - r-'k-\ (80) 

At the wavelength ^ = A '/^/o> the amplitudes become of order of the in- 
terline separation bj^ ^ Iq and the notion of bundles looses meaning — the cascade 
of bundles is cut off. The scenario at the scales A < Ab is rather peculiar and the 
next subsection is devoted to its description. 

5.3 Reconnections of adjacent lines 

In the regime of bundle crossings, in view of Eq. ( 1801 ), the Kelvin-wave am- 
plitudes are steeply decreasing with the wavenumber, so that at the wavelength 
Ab = A'/'^Zq- where the amplitudes become of order of the interline spacing Iq, 
the vortex lines are only slightly bent, bj^k <C 1 . This poses an interesting question 
of what is driving the cascade at the wavelengths below Ab- Since the amplitude 
spectrum can not change discontinuously, at the scales adjacent to Ab we must 
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Still have bf^k ^ 1, so that the mechanism of self-reconnections is strongly sup- 
pressed. On the other hand, the kinetic times of the purely non-linear regime are 
still too long to carry the flux e . We thus conclude that in some range of length 
scales Ac < A < Afo there should take place a build up of the relative Kelvin-wave 
amplitude supported by the energy flux from the scale Ab until bi^k ^ 1 and the 
self-reconnections can take over the cascade, the amplitudes bjt being defined by 
the condition of constant energy flux per unit length and the crossover scale A: 
being associated with the condition bj^^i^i^ ^ Ac- The observation crucial for un- 
derstanding the particular mechanism of this regime and thus finding bji is that 
each nearest-neighbor reconnection (happening at the rate oc A /A^ per each line 
element of the length ^ A^,) performs a sort of parallel processing of the en- 
ergy distribution for each of the wavelength scales in the range [Ac , A^] . For the 
given wavelength scale A ^ fc^', the energy transferred by a single collision is 
oc A{bk}if'k, and with the above estimate of the collision rate per the length A^, 
this readily yields the spectrum 

bk ^ loihk)-'!^. (81) 

and, correspondingly, Ac ^ l^/A^I'^. Note that the rise of the relative amplitude 
bj^k oc ^1/2 implied by Eq. dSlll describes the process of fractalization of the vortex 
lines with its culmination around the scale Ac- Thus, in this regime, the energy is 
actually contained at the low-end of the inertial range. We refer to this circum- 
stance as a deposition of energy and discuss the relevance of this term in Sec.|7] 
Although the deposition of energy makes this regime formally different from a 
standard cascade setup, its presence has no effect on the rest of the inertial range 
and thus on the overall cascade efficiency. 



5.4 To the Kelvin-wave cascade 

The stages of vortex-bundle reconnections and the reconnections between adjacent 
vortex lines essentially complete the transformation of the quasi-classical tangle 
that is necessary to connect the classical regime to the Kelvin-wave cascades on 
individual vortex lines. At the wavelengths ^ Ac the self-reconnections on the vor- 
tex lines take over the cascade and the picture at the shorter scales is qualitatively 
not different from that of the non-structured tangles. The self -reconnection regime 
continues in the range A* <C A ^ A^ with the spectrum b^^ k^^ . Kia. sufficiently 
small wavelength A*, the strongly turbulent cascade of Kelvin waves is replaced 
by the purely non-linear cascade. The spectrum of Kelvin-wave amplitudes in the 
non-linear cascade is given by Eq. ( 1531) . The value of A* can be determined by 
matching the energy flux e with / pl^, where bj^^ijx^ ^ k~^ ^ A*. With Eq. d76] |. 
we then obtain 

A*=Zo/a1/2. (82) 

At r = Kelvin waves decay emitting phonons. This dissipation mechanism 
is negligibly weak all the way down to wavelengths of order Aph, where the rate 
of energy radiation by the vortex lines becomes comparable to e. As we show in 
Sec.|6]this condition with the estimate for the power of sound radiation by kelvons 
yields 

Aph^A2V3i[»c/c/o]''/"/o. (83) 
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The scale Apj, ^ A* gives the lower dissipative cutoff of the Kelvin-wave cascade. 

The overall decay scenario is summarized in Fig. [3] where we plot the resulting 
Kelvin-wave spectrum in the whole quantized regime. 
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Fig. 3 Spectrum of Kelvin waves in the quantized regime of quasi-classical tangles. The in- 
ertial range consists of a chain of cascades driven by different mechanisms: (1) reconnections 
of vortex-line bundles, (2) reconnections between nearest-neighbor vortex lines in a bundle, (3) 
self-reconnections on single vortex lines, (4) non-linear dynamics of single vortex lines without 
reconnections. 



5.5 Scanning by finite temperature 

In this subsection, on the basis of the presented decay scenario at T = 0, we de- 
scribe a theory of the low-temperature dissipative cutoff of the cascade ITsh . We 
shall demonstrate that the temperature dependence of the total vortex-line length 
density L comes from the dependence of the wavelength scale Acuttoff > which the 
cascade ceases due to the mutual friction of vortex lines with the normal compo- 
nent, on the dimensionless friction coefficient (T 0). The main result 
of the theory is the prediction for the function InL(lna), which turns out to have 
a very characteristic form reflecting the four qualitatively distinct wavenumber re- 
gions of the cascade in the quantized regime. The results of the theory are in an 
excellent agreement with the recent experiments by Walmsley et al. Ill9ll24ll giving 
a strong evidence for the highly nontrivial scenario of low-temperature turbulence 
decay. 
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If the vortex lines were smooth, the vortex-line density L would be simply 
related to the interline separation as L = l^^. However, the presence of a fine wave 
structure on the lines can make the total length many times as large as Zq ^ (or even 
infinite in the limit of fractal lines). The increase of L due to the presence of KWs 
is related to their spectrum by IH] (cf. Eq. ( |59] |) 

/■^'cutorr{«) , 

ln[L(a)/Lo] = y_ {bkkfdk/k. (84) 

Here ^cutoff ^ 1 /■^cutoff > k is the smallest wavenumber of the KW cascade (not to 
be confused with the smallest wavenumber of the Kolmogorov cascade) at which 
the concept of a definite cutoff scale is meaningful, and Lq is the "background" 
line density corresponding to ^cutoff ^ k. There is an ambiguity in the definition of 

associated with the choice of the spectral width of the scale k, which is fixed in 
Eq. ( 184b by setting the proportionality constant on each side to unity. 

At r = 0, the cascade is cut off by the radiation of sound (at least in '^He) 
at the length scale Acutoff = ■^h given by Eq. ( [83] ). Changing the temperature one 
controls Acutoff > ■^h in Eq. (184b . which allows one to scan the KW cascade ob- 
serving qualitative changes in L{T) as Acutoff traverses different cascade regimes. 
The existence of a well-defined cutoff is due to the fact that the cascade is sup- 
ported by rare kinetic events in the sense that the collision time Tcou = Tcoii(£,^) 
is much larger than the KW oscillation period, Tper = Tper(^)- The dissipative time 
Tdis = T(iis(a,fc) ^ Tper/0£, as we show below, is the typical time of the frictional 
decay of a KW at the scale k. Thus, the cutoff condition is Tdis(«,fc) ^ Tcoii(e,^), 
which implies that the energy dissipation rate at a given wavenumber scale be- 
comes comparable to the energy being transferred to higher wavenumber scales 
per unit time by the cascade. It is this condition that defines the cutoff wavenumber 
^cutoff = ^cutoff (£ J (^)- Decreasing T and thus (X{T), one gradually increases ^cutoff > 
thereby scanning the cascade. In view of Eq. ( |84| l. this in principle allows one to 
extract the KW spectrum. 

At finite T, dissipative dynamics of a vortex line element are described by the 
equation (omitting the third term in the r.h.s., which is irrelevant for dissipation) 

mm 

s = v(s) + as' X [vn(s) - v(s)] . (85) 

Here v(r) is the superfluid velocity field, Vn(r) is the normal velocity field, s = 
s(^ , f ) is the time-evolving radius-vector of the vortex line element parameterized 
by the arc length, the dot and the prime denote differentiation with respect to time 
and the arc length, respectively. 

At a ^ 1 the superfluid and normal components are strongly coupled and the 
KWs are suppressed. In this case, the cascade must cease before it enters the 
quantized regime, i.e. Acuttoff }^ fQ. If, however, the mutual friction is small, then 
a^' ^ 1 gives the characteristic number of KW oscillations required for the wave 
to decay. Indeed, to the first approximation in 1/A, v(s) in Eq. ( 185b is given by 
the local induction approximation dl lb . In the range A ^ rg, which we will be 
interested in, the normal component is already laminar and thus the field Vn has no 
structure at these small scales. In addition, as long as a <C 1, the disturbance of the 
normal component caused by vortex line motion as these scales can be neglected. 
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Therefore, Vn(s) can be treated as a constant in Eq. ( 1851) , which makes it irrele- 
vant for KW dissipation. In this case, Eqs. dSSl l, (fTTT) give the rate at which the 
amplitude of a KW with the wavenumber k decays due to the mutual friction, 

bk^-acokbk, (86) 

while the KW dispersion is (Oi^ = {K/An)Ak^. Here and below we omit factors 
of order unity, which are subject to the definition of the spectral width of the 
wave. These factors can not be found within our theory, but can be extracted from 
experimental data (as we do it below) and from numerical simulations. Since the 
energy per unit line length associated with the wave is Ej^ ^ Kpco^bf, the power 
dissipated (per unit line length) at the scale ^ is given by 

n{k) ^ UKpa^bl (87) 

where p is the fluid density. In the following, we analyze ST as (X{T) scans 
through the regimes summarized in Fig. [3] 

Regime (1). — At this stage the vortex lines are organized in bundles; the am- 
plitudes of waves on these lines are given by Eq. (180b . One can estimate the 
total power lost due to the friction of these bundles at the wavenumber scale 
r^^ <C ^ <C ^ ^ per unit mass of the fluid as 

edis {k) = ^^n{k) aK^A/l^. (88) 

Here, the first factor in the r.h.s. is associated with the correlation volume at this 
scale ^ ' and the second one stands for the number of vortex lines in the 
volume. Note that the dissipated power is constant at all the length scales within 
this regime and, since the interline separation is related to the Kolmogorov flux by 
Eq. ( 176b . it is simply given by ae. Thus, when a ^ 1 the kinetic channel in the 
whole regime (1) becomes efficient and the cascade reaches the scale Ab, where 
the notion of bundles becomes meaningless. That is, the regime (1), as opposed 
to the regimes (2)-(4), is not actually scanned by a — its inertial range develops 
as a whole while a evolves from a ^ 1 to a ^ 1 . [In the purely theoretical limit 
of exponentially large A, when the inertial range of the regime (1) occupies many 
decades, a (loose) dissipative cutoff appears in the regime (1) as well; the cutoff 
wavelength in this case can be roughly estimated as ln(/'o/Acutoff) ^ ^/OC-] 

Regime (2). — In the range k\, <^ k kc ^ , the spectrum of KWs is given 
by bk ^ lo{k\,/ky^^ (Eq. dSOl l) which for the dissipated power yields 

edis(fc) - n{k) ^ aK^Ay^k'/h . (89) 
P'o 

The condition e^is (k) ^ £ gives the cutoff wavenumber 

^cutoff = ^lA-y^a-y^ {In/h) , fcb < fccutoff < k, , (90) 
where ^2 is some constant of order unity. Then, from Eq. ( |84b we obtain 

^ c2(;fcb/o)'[fccutoff/fcb - 1] . (91) 
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Here, we set bj^ = Clo{k\,/k)^^^, where C is a constant of order unity. The overall 
magnitude of bj( in the other regimes follows then from the continuity. 

Regime (3). — In this regime, supported by self-reconnections, the spectrum 
is given by bk ^ (up to a logarithmic prefactor). The corresponding energy 
balance condition yields the cutoff scale: 

^cutoff = 5 {Aay^^^ {2n/lo), K < ^cutoff < K ■ (92) 

With the logarithmic prefactor taken into account, Eq. ( l66l) , the spectrum in this 
regime reads b]^ = C[\ +c\ ln(^/fcc)]~'^^(V^c^b/^)fe' where C3 is a constant of 
order unity. Then the relative increase of vortex line density through this regime 
is given by 

^(fccutoff) 

L{h) 

where v = C'^k^l / Cy 

Regime (4). — Since the spectrum of the purely nonlinear regime, bk ^ k^^l^, 
(Eq. ( 153b ) is steeper than the marginal b]^ ^ k^^ meaning that the integral in 
Eq. ( 184b builds up at the lower limit, as soon as ^cutoff ^ K- the line density 
i(^cutoff) starts to saturate and becomes independent of ^cutoff at ^cutoff > K- The 
energy balance gives the dependence fccutoff(o')i 

^cutoff = ^4 A a-^l^ (27r/Zo) , ^cutoff > h , (94) 

where (^4 is an unknown constant. The coefficient in the KW spectrum is fixed by 
continuity with the previous regime, yielding 

bk = C[\ +cl\n{k^/k^)]-''l^-^/k^ky\-^IHo. Eq. IH thus yields 



1 + cUn 



^cutoff 



(93) 



^(Wf) {5c^mMl 

L{k,) ~ l+c2ln(^*Ac) 



2/5 



\ ^cutoff / 



(95) 



The continuity of ^cutoff leads to the following constraints on the free coeffi- 
cients in Eqs. (go}, and dH, ^3 = & A^l^al'^, ^4 = ^2 A^l^ O^'^ al'^ , 
where A is the value of the friction coefficient at which the cascade is cut 
off at the crossover between the regimes (2) and (3) and a* ^ 1/A^ corresponds 
to the one between (3) and (4). Introducing ab<\, corresponding to the crossover 
from the regime (1) to (2), we rewrite Eqs. ( l93] l, and (|95] l in terms of a: 



In — ^ = A2 



(ab/a)'/3-l 



L{a) _ 
L(a] 



l + {ci/2) In^j 



L(a*) 



I+A4 



1-1 « 

Of* 



1/41 

, a<ca*, 



(96) 
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where 

V = {2nf^id/A'/^ia,ai,y/'cl (97) 




c b 

Fig. 4 (Color online.) Fit of the experimental data (closed circles and triangles) adapted from 
Ref. l24ll (/o » 5 X IQ^^cm giving A w 13). The high temperature measurements of Refs. IsTll . 
fish are represente d by open squares and open triangles respectively. The form of a{T) is taken 
according to Ref. i44[| at J > 0.5K (roton scattering) and according to Ref. l35ll at J < Q.5K 
(phonon scattering). The fitting parameters are Aj ~ 0.25, A4 0.26, C3 Ri 3.0, v ~ 0.14 with 
Kb ~ 3.5 X 10-2, ^ 2.5 X lO"'', a, ~ 2.0 x 10"^. 



Although the excellent quantitative agreement with the experiments can be 
sceptically attributed to the large number of fitting parameters, it is remarkable 
that the nontrivial qualitative behavior of L{a) exhibits the form peculiar to our 
decay scenario. As a decreases from a 1 to the values significantly smaller 
than unity, the line density L increases only by some factor close to unity 1.5 
in the experiment), which reflects the formation of the regime (1) driven by the 
reconnections of vortex bundles. During the crossover from the region (1) to (2), 
the increase of L is minimal — a shoulder in the curve L{a) arises. It is only well 
inside the region (2) that the increase of L becomes progressively pronounced, 
and, at the crossover to the region (3), the function L{a) achieves its maximal 
slope, determined by the developed fractalization of the vortex lines necessary to 
support the cascade within the interval (3). As the cutoff moves along the interval 
(3) towards higher wavenumbers, the slope of L{a) becomes less steep due to 
the logarithmic decrease of the characteristic amplitude of KW turbulence. When 
the cutoff passes the crossover to the regime (4), the curve L{a) gradually levels. 
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0.01 r 




Fig. 5 Spectrum of Kelvin waves in the quantized regime quantified, apart from the regime (1), 
by the fit to experimental data, Fig.|4] In view of Eq. <84t , the constants C and <^2 can be found 
only in the combination (^oC ~ 0.049. The regimes (l)-(4) correspond to those in Fig.[3] 



[Note that because of the closeness of the KW spectra in the regimes (3) and 
(4) the crossover between them, as expected, is not distinct, which results in a 
significant uncertainty in the choice of a*.] Fitting the experimental data fixes the 
values of the dimensionless parameters, thus revealing the quantitative form of 
the KW spectrum shown in Fig. |5] Note, however, that in view of the fact that the 
experimental A ^ 10 is not so large we would expect certain systematic deviations, 
but a significant scatter of the experimental data does not allow us to assess them. 



6 Kelvon-Phonon interaction 

6.1 Hydrodynamic action. Hamiltonian of Kelvon-Phonon interaction 

The general problem of kelvon interaction with other modes is relevant far beyond 
the scope of superfluid turbulence. Since the early study of phonon scattering by 
vortex lines ll45ll . much interest has been attracted by the problems of the inter- 
action of Kelvin waves with density-distortion modes H [S EE El, IH SI]. In 
neutron stars, the excitation of Kelvin waves due to the interaction with the nuclei 
in the solid crust is suggested to be the main mechanism of pulsar glitches ll4lll . 
Nonlinear kelvon dynamics is also becoming an attractive topic in the field of 
ultra-cold gases, where one can study kelvon coupling to the modes of the Bose- 
Einstein condensate in situ, as, e.g., in the problem of kelvon excitation by the 
quadrupole mode lElElEl] 

In this section we describe a systematic approach (developed by the authors in 
Ref. iflSll ') to the problem of interaction of phonons with vortices in the hydrody- 
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namic regime, i.e. when any physical length scale is much larger than the vortex 
core size ao, which allows one to describe vortices as geometrical lines ijH]. We 
derive the interaction Hamiltonian basing the analysis on the small parameter 

j3=flofc<l, (98) 

where k is the largest wave number among kelvons and phonons. To employ the 
transparent description in terms of the normal modes, we confine ourselves to 
the case of weak nonlinearity, which for kelvons implies the small parameter a, 
Eq. (|30l l, meaning that the amplitudes bf, of the Kelvin waves of the typical wave- 
length ?i ^ are much smaller than ?i, ai(= bj^k <C 1. For phonons this requires 
that 77 <C M, where T] is the number density fluctuation in a sound wave and n is 
the average number density. The obtained result allows us to rigorously describe 
the radiation of sound by kelvons, which we apply to the problem of superfluid 
turbulence decay at zero temperature in the next subsection. 

Let us first describe the qualitative ideas behind the derivation. The condi- 
tion ( |98l l implies that the typical vortex-line velocities are much smaller than the 
speed of sound. Along with 77 ^ « it leads to the fact that the vortex-phonon cou- 
pling contributes only small corrections to the dynamics of the non-interacting 
vortex and phonon subsystems. Therefore, a perturbative approach is applica- 
ble, provided the interaction energy is written in terms of the canonical vari- 
ables. Normally, the form of the canonical variables comes from the solution 
of the interaction-free dynamics. However, when studying vortices separate from 
phonons, one naturally neglects the compressibility of the fluid, which we did in 
Sec.|2] since finite compressibility leads only to higher-order "relativistic" correc- 
tions to vortex dynamics. As a result, the dynamics of vortices are described by 
the Hamiltonian ( 129b , written in terms of the geometrical configuration of the vor- 
tex lines. When the vortices are absent, the phonon modes come from the bilinear 
Hamiltonian for the density fluctuation 77(r,f) and the phase field (p{r,t), which 
determines the velocity in the density wave (see, e.g., lli^ ). If finite compress- 
ibility of a superfluid is taken into account in order to join the subsystems, the 
positions of vortices and the fields T] {r,t),(p{r,t) are no longer the sets of canon- 
ical variables because of the variable-mixing term in the Lagrangian. This makes 
this perturbative problem quite peculiar since introducing the interaction, one has 
to simultaneously reconsider the canonical variables. 

One can easily see that the standard vortex parametrization used in Sec.|2]fails 
to capture the physics of vortices when phonons are present. In Eq. (|4]l, the geo- 
metrical position of the vortex lines unambiguously determines the instantaneous 
velocity field configuration at arbitrary distances, which is inconsistent with the 
finite velocity of propagation of excitations in a compressible medium. From this 
simple physical argument, we can already guess the form of the canonical vari- 
ables: the long-distance part of the velocity field produced by vortices in Eq. (|4| 
should actually belong to phonons. 

The small parameters allow us to obtain an asymptotic expansion of the canon- 
ical variables by means of a systematic iterative procedure. Physically, the proce- 
dure restores the retardation in the adjustment of the superfluid velocity field to the 
evolving vortex configuration. It qualitatively changes the structure of the Hamil- 
tonian with respect to the terms responsible for the radiation of sound and the 
relativistic corrections to the vortex dynamics. 
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Hydrodynamic Lagrangian: standard parametrization. Long-wave superfluid 
dynamics at zero temperature are described by the Popov's hydrodynamic action 



S= dtd^r 



-{n + ri)4> 



2mn 



2>f 



(99) 



Here the spatial integral is taken over the macroscopic fluid volume, 'P{r,t) is the 
phase field, which determines the velocity according to v = {l/mo)'V<P (h = 1), 
niQ is the mass of a particle, x is the compressibility, the dot denotes the derivative 
with respect to time, and the microscopic length scale is defined by aq = \/ >c/nmo. 
The vortex core radius is of order ao- 

The phase is non-single-valued and contains topological defects, the vortex 
lines. The velocity circulation around each vortex line is quantized; 



<j)y^{r)-dr = 2n. 
The defects can be separated from the regular contribution: 



(100) 



(101) 



where 4>o is non-single-valued and satisfies ( 1 100b . while 9 is regular and V(p is 
circulation-free. The standard decomposition into vortices and phonons is done by 
introducing an additional constraint that to the zeroth approximation eliminates 
the coupling between 4>o and 9 in the Hamiltonian, namely 



0. 



(102) 



(Physically, this parametrization is suggested by the velocity potential of an in- 
compressible fluid.) With Eqs. ( |101| ),( fTU2l ) the Lagrangian becomes 



Jd^r[-n^Q-ri<p-ri4>o\ -H, 



where H = //„ 



H' 

^^int' 



2mo J 



ffph = jd^ 



" 2mn 



d^r 



77|V4>o| +277V9-V^>o 



(103) 

(104) 
(105) 
(106) 



The coupling between the vortex variable, <i>Q, and the density waves, {77,9}, is 
determined by H[^^ and the time derivative term Jd^rri<pQ, both being first-order 
corrections to the non-interacting parts. 

Noninteracting case. Following standard perturbative procedure, we first ne- 
glect the coupling between the vortex and phonons to find the non-interacting nor- 
mal modes. The vortex part of the Lagrangian is then given by Lyor = ~n Jd^r (pQ — 
//vor- For the sake of simplicity, from now on we consider a solitary vortex line; 
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the generalization is straightforward. Let the two-dimensional vector Po(zo) = 
{xo{zo) ,yoizo) ,0) describe the position of the vortex line in the plane z = zo of a 
Cartesian coordinate system, where the z-direction is chosen along the vortex line. 
The field is a functional of Pq{z), hence 

Jd^r <PQ = jdz po • ^ jd'^r <Po. (107) 

To obtain / d^r 4>q, it is sufficient to calculate the integral / d^r 5^q, where 5^q 
is the variation of the phase field due to the distortion of the vortex line by 5Po(z). 
Using the identity 4 (p^) = 4, with p = (jc, v,0), and Eq. ( 1102b . obtain 

jd^r84>o=^-jd^ry-[84>oy{p^)]. (108) 

The variation 5'S'o(r) can be viewed as being produced by two vortex lines with 
opposite circulation quanta and separated by 5po(z). In view of Eq. ( 1 1001) the field 
54'o(r) experiences a jump of 2n across the surface that extends between these 
vortex lines along the vector field 5pg(z), therefore the integration volume must 
have a cut along ,5^ . Applying the Gauss theorem, we rewrite Eq. ( 11081) as the 
surface integral (1/2) §y 5^o (p • dS) yielding 

Jd'r <Po = n jdz [i X Po(z)] •Po(z). (109) 

Introducing the complex variable w(z) = y/nmoK/2 [x{z) + (v(z)], where K = 
2n /mo is the velocity circulation quantum, obtain 

ivor = ^ Jdz [iw*W~ iw* w] - Hyor [w,w*], (110) 

which implies that w(z) and w*{z) are the canonical variables with respect to Lyor- 
The energy ( |104| l. rewritten in terms of w(z) and w*{z) gives the vortex Hamil- 
tonian jlOi . or the pseudo-Hamiltonian ( 1291 ) with an appropriate regularization. 
We need only the bilinear term of the expanded with respect to a*; ^ 1 Hamilto- 
nian since the kelvon interactions are irrelevant for our problem. We change the 
notation for the kelvon spectrum reserving (0^ for phonons: 

Hyoi « J^Skala,, , = ^ [ln{l/kao)+CQ]k^, (111) 

k 

with and aj, defined in Sec.|2] The sound waves are described by the Lagrangian 
Lph = / d^r{-r\<p) -//ph[T], 9], with //ph given by dlOSb . We assume that the sys- 
tem is contained in a cylinder of radius R with the symmetry axis along the z- 
direction and that the system is periodic along z with the period L. In the cylindri- 
cal geometry, the phonon fields T](r, 0,z), (p{r,d,z) are parametrized by phonon 
creation and annihilation operators c,, c] as 



(p = -/£ ^/lJ2^ [XsCs- X^c] ] , (112) 
Xs = Xs{r,9,z) = J%,ng,{r) Y,„{e) ^,,(z), 
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where ^„„/,(r) = {nqr/RY'^JnMrr), Y,n{e) = (27r)-i/2exp(/m0), J^^^(z) =L-^I^ cx^[iq^z), 
s stands for {qr,m,q^}, and Jm{x) are the Bessel functions of the first kind. The 
phonon Hamiltonian then reads 

^^ph = (Osdc,, (Os = cq, (113) 



with q = y^q^'+q^ and c = y^n/xmo = 1 /moflQ. 

The effect of coupling: change of canonical variables. Now we address the 
coupling between phonons and vortices. In terms of the obtained variables, the 
Lagrangian ( 1103b takes on the form 

L = Y,iakal+Y,iCscl~T-H, (114) 

k s 

where T = Jd^r 77<j>o = (^s,Cs,c],cl}, and H = H ^p, + H^h + H^^^^. 

The coupling term T plays a special role in the Lagrangian (II 14b . This term is lin- 
ear in time derivatives of the variables and thus can not contribute to the energy in 
accordance with the Lagrangian formalism. Moreover, because of the time deriva- 
tives in T, the equations of motion in terms of {ajc,al},{cs,cl} take on a non- 
Hamiltonian form. This implies that the chosen variables become non-canonical 
in the presence of the interaction, and therefore the total energy, H, in these vari- 
ables can not be identified with the Hamiltonian. 

There must exist such a variable transformation {aj^ , } , {c, , c] } — > {a^ , a], } , {c, , c] } 

that restores the canonical form of the Lagrangian, L = ^.k^^k^l +L.s''?.sc] — 
H{af;,a^i^,Cs,cl}. The canonical variables are obtained by the following iterative 
procedure. The term T is expanded with respect to 1, j3 ^ 1 and 7] <C « 

yielding T = T^^^ -f T*-^^ H . Then the variables are adjusted by aj^ — > % + a[^', 

Cs Cj + Cj^', where fl['^({a<:,aj,Cs,Cj}) and ci^\{ak,al,Cs,cl}) are chosen to 

eliminate the term in ( 11141) . As a resuk, T ^ + T'^^' H , where the 

prime means that the structure of the remaining terms has changed. At the next 

step, r''^' is eliminated by a^. aji + a^^\ Cj and so on. By con- 
struction, the canonical variables are given by 5^ = a^. + +"^1^^ > = 

Cs + Cj^' + cp' H , and likewise for the conjugates. In practice, only the first few 

terms are enough, as the rest ones give just higher-order corrections. 

The explicit expression for T is obtained following the steps of the derivation 
of Eq. ( 1109b . The only difference here is that the role of the auxiliary function 
is played by Q, defined by A Q{r) = 7] (r) : 

T = 2n ldz[ixVQ{Po{z),z)]-i3oiz). (115) 

In accordance with qrpo ^ jSfcpo ^ 1, we expand the radial functions, retaining 
only the greatest term for each particular angular momentum m. Switching to the 
cylindrical coordinates, Pq = (po cosy, po siny, 0), and noticing that 



rPo 3- - — 3^ 
dr Po dd 



r=po,e=7 

d{wn/dt, m>0 

d{w*\'"\)/dt,m<0 ' ^ ' 
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obtain 



-iA 



-iBs.k\...k„,Cs ■Jji'^k, 
where the sum is over all s with m j^O and 



H.c, 



(117) 



B 



s.k. 



= (-l)H0(-,«-l)A,5; 



\/q/moi 



2l'"l/2+i |m|! 

where ©(m) = < n' ^ n ■ Thus, 
^ ' [ 0, m < 



L(l'«|-l)/2/?l/2 



^-2 



(118) 



(119) 



(we omitted the terms that do not contain phonon operators and thus result only in 
relativistic corrections to the kelvon spectrum and kelvon-kelvon interactions), 



c, = c. 



(1-Vo) I 

ki...k„^ 

+Bs.ki...k,„aki ---ak,,, 



, s = {qr,m,q^}. 



(120) 



The interaction Hamiltonian. The Hamiltonian H is given by the energy ( 1104b - 
(IIO6I 1 in terms of the variables {fl|t,a^},{c5,cj}. Up to neglected relativistic cor- 
rections, the variable transformation does not change the spectrum of the elemen- 
tary modes: the zero-order Hamiltonians are given by (II 1 lb . dl 13b in terms of 
{c^jc]}. The transform ( 1120b applied to ( 1113b generates the interaction 

term 



(rad) 



^1(1 

s,{ki} 



co.A 



'sAsl|.../t, 



k\ 



(0,B 



s,kt ...k,„ <2ii 



■^kjl 



■^k 



H.C. 



(121) 



Remarkably, the energy term oc J t/^r?] | V<J>o|^ in (1106b . which results in the same 
as ( 11211) operator structure, is irrelevant, being smaller in /3 <C 1 . It can be checked 
straightforwardly, that the term oc Jd^r r]V(p ■ V^Pq in Eq. ( 11061) gives Fetter's 
amplitudes of the elastic and inelastic scattering of phonons [4511 (see, however. 
Sec. 16.31 ). In addition, this term leads to a macroscopically small splitting of the 
phonon spectrum due to the superimposed fluid circulation. 

A kelvon carries a quantum of (negative) angular momentum projection. The 
interaction ( 1121b explicitly conserves the angular momentum: a real process of the 
emission of a phonon with the angular momentum (— m) requires an annihilation 
of m kelvons. 
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6.2 Sound emission by Kelvin-wave cascade 

Since e^^ ^ {aok)(iOj^, the total momentum transferred to phonons in a radiation 
event should be small in order to satisfy the energy conservation. Thus, the radi- 
ation by one kelvon on an infinite line is kinematically suppressed. The leading 
radiation process is the emission of the m = —2 (quadrupole) phonon mode, the 
events involving more than two kelvons being suppressed by ^ 1 . First-order 
processes of the two-phonon emission come from the term oc J (fir T] V(p • V4>o of 
(1106b . The amplitude of these processes is suppressed by the relativistic parameter 
j3< 1. 

The qualitative observation that the sound emission is due to the quadrupole 
radiation is already sufficient to restore the formula for the power iT^, radiated by 
the Kelvin waves at the wavenumber scale k per unit vortex-line length. From gen- 
eral hydrodynamics [fsTt], njt must be proportional to the square of the third-order 
time derivative of the quadrupole moment, oc e^b\b^_j^, and inversely proportional 
to the fifth power of the sound velocity c. The rest of the dimensional coefficients 
are restored unambiguously, since the only remaining time and length scales are 
set by K and k and the dimensions of mass come from the fluid density p = nmo, 
which yields 

O £^k 

This result is valid up to dimensionless prefactors of order unity. 

We can use the Hamiltonian ( 11211) to obtain an accurate formula for the turbu- 
lence decay rate due to the sound radiation. The kelvon occupation number decay 
rate is given by hj^ = — '£,s.ki ^s.kMi^ where Ws.k.k^ is the probability of the event 
\0s, rijc, n|c^) \ ls, n^ — l, nic^ — l) per unit time. Applying the Fermi Golden Rule 
to Ws^k,ki with the interaction (1121b and replacing the sums by integrals, obtain 

= - ^-^j^ [ln( 1 /a,k) + Co] ' {k/cf k' nj. (123) 

The total power radiated by kelvons per unit vortex-line length is then obtained 
from J7t = — ^k'~k ^k''^k' which allows to restore the dimensionless prefactors 
in Eq. ( 1122b . However, since the theories of non-structured and quasi-classical 
tangles can not trace coefficients of order unity we shall confine ourselves to using 
Eq. (1122b for determining the cascade cutoff scale Apf,. 

In both non-structured and quasi-classical tangles, at the high wavenumbers 
where sound radiation becomes appreciable, turbulence decay is due to the purely 
nonlinear Kelvin-wave cascade with the spectrum given by Eq. (153b . The condition 
for finding the cutoff scale Aph is Tlj^^^^ ^ 6. Clearly, the value of 9 depends on the 
structure of the tangle. In non-structured tangles 0ns ^^P^^/^o^ where Iq is the 
interline separation and A = ln(Zo/«*), which gives 

Aph ^ A^^l^^ [K/clof^l^^ lo (non - structured tangles). (124) 

For the quasi-classical tangles 0qc ^ £plo, which with the Eq. 1761 gives dqc 
K^PA/Iq ^ 6„s/A. Thus, for the cascade cutoff we obtain 



Aph - A^^/^' [jc/c/o]^^''^' lo (quasi - classical tangles). (125) 
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A comment is in order here. In superfluid turbulence, kelvons do not live on 
infinite vortex lines, as we assumed throughout this section, but are superimposed 
on vortex kinks with a typical curvature radii Rq much larger than the kelvon 
wavelength, Rq^ k~^. This circumstance led Vinen [10] to obtain a dimensional 
estimate 17^°=^^°= «t for the power radiated per unit length of the vortex line, 
which is qualitatively different from our Eq. (11221) . In the quasi-particle language, 
n'f, implies that the radiation is governed by the conversion of one kelvon into a 
phonon. Vinen argues that this process becomes allowed in due to the finite size 
of the kinks, or, equivalently, the kelvon coupling to a kink lifts the ban on single- 
kelvon radiative processes by effectively removing the momentum conservation 
constraint. We note that the probabilities of such elementary events are likely to 
be suppressed exponentially, as it is generically the case, say, for soliton-phonon 
interactions (see, e.g., Ref. ll52h and references therein.) The processes involving 
kelvon-kink coupling should contain an exponentially small factor ^ exp(— 7?ofc), 
which arises from the convolution of the smooth kink profile with the oscillating 
kelvon mode. 



6.3 Elastic and inelastic phonon scattering 

The main advantage of the Hamiltonian formalism developed here is that it allows 
to reduce dynamics to elementary scattering events tractable by the standard scat- 
tering theory, thereby providing a general systematic and physically transparent 
theoretical tool. So far, the Hamiltonian framework allowed us to address the most 
important remaining problems of vortex dynamics — reconnection-free kinetics of 
Kelvin waves and phonon radiation by Kelvin waves. The purpose of this section 
is to use the developed formalism to revisit the longstanding problem of phonon 
scattering on vortex lines. The interest to this problem has been maintained by 
the controversy on the value of the elastic scattering cross-section and the corre- 
sponding force between a vortex and the normal component (see Ref. IsbIi and 
references therein). The elastic scattering cross-section was originally obtained by 
Pitaevskii [54] and later confirmed by Sonin (Ref. |53] and references therein) bas- 
ing on hydrodynamic equations of motion. However, the Hamiltonian approach to 
this problem suggested by Fetter ll45ll (and rederived later in Ref. ifSSll ^ gave a dif- 
ferent result. Within the formalism developed in this section it is straightforward 
to see, as we do here by means of a direct calculation, that the discrepancy is due 
to the non-canonicity of the standard vortex and phonon parametrization used in 
the Hamiltonian of Ref. ll45h . The result of Ref. ll45ll actually corresponds to the 
case of a pinned vortex, which was first noted by Sonin IsbIi . and the effect of vor- 
tex motion during the scattering is properly accounted for by the transformation 
to the canonical variables. 

We also discuss inelastic phonon scattering and show that the only available 
to our knowledge solution of Ref. 145] suffers from the same problem: the non- 
canonicity of the variables overlooked there makes a non-trivial contribution to the 
scattering cross-section. The accurate result can be obtained from a finite number 
of diagrams presented here (Fig. [H, but in view of the length of the resulting 
expression and absence of its immediate application, we do not derive it in a closed 
form here. 
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Vortex-phonon Hamiltonian in the plain-wave basis for phonon fields. For 
our purposes, it will be convenient to rewrite the vortex-phonon Hamiltonian of 
Sec. 16.11 in the plane-wave basis for the phonon fields explicitly obtaining the 
terms responsible for phonon scattering: 

77(r)=£.y^W2V + 4] , 

q 

•PW = v/l/2Vo;q»c[e''"-Cq-e-'Vct] , (126) 
q 

where a)^=cq and V = is the system volume. The (harmonic part of) phonon 
Hamiltonian is given by Eq. (11131 ) with s labeling different wavenumbers q. 

A calculation analogous to that of Sec. l6.1l vields the relation for the canonical 
variables in the form 

ak = dk + Ykki q [cq + cLq] + smaller terms, 

/ \ 1 /2 ^2 ^ ^2 

Ykki q = - ^77 2 ' ^ki+q„k, (127) 

(we omitted the terms that do not contain phonon operators and thus result only in 
"relativistic" corrections to the kelvon spectrum and kelvon-kelvon interactions) 
and for phonon operators 

Cq = Cq + a*,, dk + Gtik + '/iqA-i ^/^i d^, + i/jq^, ^2"-^ ^-k2 ' + smaller terms 



v,^ / \ 1/2 + ''?v c- 



l^qk\k2 



(yqk 

' (32«y)i/v 



^ki+k2,qz 



(128) 



Here and below we assume summation of products over repeating indices. 

As in the cylindrical case, up to neglected relativistic corrections, the vari- 
able transformations ( |127| l. ( |128l l applied to Eqs. dl 1 lb . ( |113l l do not change the 
spectrum of the elementary modes, but generate a non-trivial contribution to the 
interaction Hamiltonian in higher orders: 



(1) 



where 



^qk ^q ^k + SqkCqdf^ 



iM^kik2^l^kidk, + iM^kik2c\d'^_,,^d'^_k_^ 



Gkki q{^k Cq + 5^ fl/ti C_q 



H.c. + smaller terms, (129) 



Sqk = (^q 5q,q' Cq'^^, 
Mqkik2 = ^q,q' Mq'/ti/tT' 
Gkkiq = £k Sk,k' Yk'kiq- 



(130) 
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vortex line 



Fig. 6 Geometry of the elastic phonon scattering. 



Another relevant contribution to the interaction Hamiltonian comes from the 
vortex-phonon coupling energy, Eq. ( 1106b (equivalent to the one suggested by Hall 
and Vinen 15&\ ). 



H 



(2) 



jdr T]V(p-Vv, 



(131) 



where Vy = 'Vv(r) is the velocity field produced by the vortex line. Taking into 
account this term alone leads to Fetter's effective interaction Hamiltonian ll45ll 
(the terms oc acc, ac^c^ are omitted) 



„{2) _ 



<ll<l2 



q2'^qi'^q2 



^<:qiq2 ^I^qi^qz 



-H.c. + smaller terms, 



(132) 



where 



Q 



R 



qiq2 



iK [W;^ (qi X q2 - z) 

1 2 "ri.^j' 



2L^\I (Oq^ Iqi-qal 



z[(q2-qi)-(x + /y)] +^(x + /y) 



qiqi 



^**1 + V^*l2j x(q2-qi) 

I 12 ■ (133) 



The complete interaction Hamiltonian up to the highest relevant order in jS <C 1, 

< 1, and |t]|/« < 1 is given by Hint = H^^i +^1^1- ^°te that in this form 
the Hamiltonian conservers the momentum along the vortex line (along the z- 
axis), but not the transverse momentum. The change of transverse momentum 
during the scattering is transferred to the vortex line as a whole (an analog of the 
Mossbauer effect), which results in a macroscopically small displacement of the 
vortex line from its initial position. The corresponding Goldstone mode can be 
straightforwardly taken into account if necessary. 
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Fig. 7 Diagrammatic representation of the transition amplitude of elastic phonon scattering on 
a vortex line. The wavy and straight lines represent phonons and kelvons respectively. 



Elastic phonon scattering. We consider the scattering geometry shown in Fig.|6] 
The vortex line has a distribution of kelvon occupation numbers {«<:}. The elastic 
scattering differential cross-section can be found as 




(134) 



where Wqj^qj is the probability per unit time of a scattering event of a phonon 
with the wavenumber qi into the one with q2, such that qi = q\ (sin0, 0, cos 0), 
<J2 = ^2 (sin0cos(/), sin sin (/) , cos ) . The probability Wq,^q, is given by the 
Golden rule, Wq,^q, = 27r5(ft)q, — WqJIVqj^qjP- To the first (Born) approxi- 
mation the transition amplitude V^^^^ is composed of the elementary processes 

shown diagrammatically in Fig.|7] Note that, in Ref. ll45ll only the term Vq^l+q^ was 
taken into account. The remaining terms describe a scattering process in which the 
vortex line undergoes a virtual deformation (produced by a kelvon with k = qi^). 

It is easily seen that physically Vql' account for the oscillatory motion of the 
vortex line induced during the scattering: if the vortex were pinned by an external 
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potential these ternis would have to vanish. The transition amplitudes in Fig.|7]are 
given by 

^qi^q2 ^^qiq2+^q2qi = 
iK ^^Ig^gi +g2)sin^esin(/) ^ 
IL"^ qj+qj- 2qiq2{sm^ 6 cos (j)+cos^ 6) ^'''^2' 

<^q. = 77^K + i)> (136) 
^qAq2= \ {n-k + i), (137) 
^qi-q2 = — 7— (138) 

Virl,. = ^^«.. (139) 

Note, that at small scattering angles the amplitude Vq"i*q, diverges and the pertur- 
bation theory is not applicable. However, the small-angle scattering does not play 
an important role for the dissipative component of the mutual friction force. Since 
£k/(Oq ^ f o^iz ^ j3 <C 1 for 1^1 = = \q2z\, we can neglect the kelvon energy in 
the denominators of Eqs. ( I135I )-( I139I ). Thus, the terms proportional to {n^} cancel 
out in the total transition amplitude. 



^qi^q2 — ^qiq2 



5'q, k 5'q2 



(140) 



meaning that alternatively the amplitude can be obtained by retaining only the 

terms of type Vq"-^q2 but evaluating them in kelvon vacuum. Substituting the ver- 
tex expressions into Eq. ( 1140b and doing the algebra yields 

iK sin^0sin(/)cos(|) + cos"0sin(/) 
2L2 1 - cos (/) 

Finally, replacing the sum in Eq. (1134b by the integral, 

da I \ [ l2dq2L?' 



'^i^q2 = TFT i T^ZTZ "qumz- \^^^) 



d^ 

we arrive at 



^ = (V/c)y^^5(0;q,-0;qJ|Vq,^, 



2 

qi^q2i ' 



(142) 



da [sin^ sin(i) cos(i) +cos2 sin(i)] ^ 

= T a - - - ('1431 

d(p ^ Snc^- (l-cos(/))2 ■ ^ 



For 9 = 7t/2 this result gives the differential cross-section obtained in Ref. 05311 . 
Due to the divergency at (/) ^ 0, this formula is only applicable in the cases where 
the small-angle scattering is irrelevant, as, e.g., in a calculation of the dissipative 
mutual friction force, which immediately follows from Eq. ( 1143b . In contrast, the 
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Fig. 8 Inelastic phonon scattering. Tiie diagrams give the total transition amplitude. The wavy 
and straight lines represent phonons and kelvons respectively. The virtual kelvons are described 
by vacuum propagators. 



calculation of the non-dissipative component is essentially non-perturbative and is 
discussed in Ref. IsbIi . 

Inelastic scattering. Employing the same formalism, we describe the inelastic 
phonon scattering, i.e. scattering accompanied by generation of Kelvin waves. In 
the diagrammatic language the transition amplitude Vqi^q,,*:, where k is the kelvon 
wavenumber, is shown in Fig. [8] Note that the extra diagrams with reversed prop- 
agators are dropped, which is accounted for by evaluating the remaining ones over 
kelvon vacuum. The scattering matrix element reported in Ref. ll45h corresponds 
to the first term oc T^qjq^ only. 



7 Conclusions and discussion 

The zero-temperature limit of superfluid turbulence is non-trivial and instructive. 
A number of small and large parameters — Eqs. ([T|, ( 130b , ( 1981 ) — and (approxi- 
mately) conserving quantities (energy, momentum, number of kelvons, LIA con- 
stants of motion), characterizing the phenomenon, bring about rich physics on one 
hand, and provide a solid basis for applicability of the analytic tools on the other. 

The central part in the theory is played by the large parameter A, Eq. ([ij. 
Implications of the condition A ^ 1 include (i) universality of the answers for 
all U(l)-type superfluids, (ii) applicability of LIA, (iii) features following from 
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(ii), of which the most important ones are (a) the integrability of LIA suppressing 
the efficiency of the pure Kelvin-wave cascade and (b) preemptive character of 
the self-induced motion of the vortices in the smallest Richardson-Kolmogorov 
eddies, responsible for cutting off the classical-fluid regime before it could expe- 
rience a potential bottleneck associated with (a). 

The energy and momentum conservation renders impossible the pure vortex 
ring cascade (originally proposed by Feynman) in which superfluid turbulence is 
supposed to decay into small rings in such a way that each ring produces smaller 
ones independently of the rest of the vortex tangle. The absence of pure Feyn- 
man's cascade and suppressed efficiency of the pure Kelvin-wave cascade create 
specific circumstances under which there exists a range in the kelvon wavenumber 
space — with the size controlled by A — where the Kelvin-wave cascade is driven 
by this or that type of vortex line reconnections. The reconnections lift the integra- 
bility constraint of the LIA and push the vortex-line length to higher and higher 
wavenumbers. 

In the theoretical limit of really large A (to be taken with a grain of salt in view 
of realistic A < 15 in ^He), we predict three distinct types of reconnection-driven 
Kelvin-wave cascades. The first two ones (in the down-the-cascade order) — (i) 
vortex-bundle-reconnection-based and (ii) neighboring-line -reconnection-based — 
apply only to the Richardson-Kolmogorov quasi-classical regime, while the third 
one, driven by local self-crossings, applies to non-structured tangles as well, in- 
cluding superfluid turbulence sustained by vibrating objects. At large-enough wave- 
numbers, the local-self-crossings-driven cascade naturally crosses over to the pure 
Kelvin-wave cascade, as soon as the power of the latter becomes sufficient to com- 
pete with the former, eventually suppressing the reconnections by rendering the 
amplitude of Kelvin waves small compared to their wavelength. 

At r = 0, the inertial range of the pure Kelvin-wave cascade is controlled 
by the small parameter j3, defined in Eq. ( |98b , that guarantees irrelevance of 
compressibility and/or kelvon-phonon interaction before the cascade reaches the 
wavelength scale A ^ Aph, where phonon emission becomes dominant, cutting 
the cascade off. In a typical case, l/jS '^ln(l/j3), implying that the inertial 
range of of the pure Kelvin-wave cascade is significantly larger than that of the 
reconnection-driven cascade(s). 

Kinetics of the pure Kelvin-wave cascade are perturbative due to the asymp- 
totically vanishing amplitude-to-wavelength ratio. This allows one to introduce 
a kinetic equation with the collision term associated with kelvon scattering. Be- 
cause of the conservation of the total number of kelvons (which is nothing but 
the conservation of the angular momentum) and one-dimensional character of the 
problem, the leading kinetic process is the elastic three-kelvon scattering. As an 
effect of the LIA integrability constraint, the kinetics are entirely due to the inter- 
actions going beyond the LIA: the local contributions to the effective three-kelvon 
scattering amplitude exactly cancel each other leaving no order-A collision terms 
and suppressing kinetic rates. 

The small amplitude of Kelvin waves at the dissipative cutoff wavenumbers 
and weak kelvon-phonon coupling guaranteed by the parameter /3 imply that the 
phonon emission is a perturbative process that can be naturally described within 
the Hamiltonian formalism. Our derivation of the kelvon-phonon Hamiltonian re- 
vealed a subtlety of non-canonicity of the prime harmonic variables with respect to 
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the interaction Hamiltonian. Adequately accounting for this fact at the level of hy- 
drodynamic action leads to a proper kelvon-phonon interaction Hamiltonian, con- 
taining terms that would be missed otherwise. With this Hamiltonian the problem 
of the phonon cutoff of the Kelvin-wave cascade becomes rather straightforward. 
The leading process is two-kelvon inelastic scattering converting two kelvons with 
almost opposite momenta into one phonon. As a by-product, we resolved some 
long-standing controversies of the kelvon-phonon interaction, showing that these 
were due to the above-mentioned lacking terms in the interaction Hamiltonian. 
[We also revealed the role of a hidden phonon-vortex interaction in the theory of 
Berezinskii-Kosterlitz-Thouless transition llssll . which is not covered here.] 

Having at hand a complete theoretical scenario for superfluid turbulence at 
r = 0, it is important to check it against experiment and numeric simulations. 
Discussing first the more simple case of non-structured turbulence, we note that a 
circumstantial evidence for the Kelvin-wave cascade — the absence of temperature 
dependence of the relaxation time superfluid turbulence in "^He at T < lOmK — 
has been observed in both experiment 10] and simulations J^. The simulation of 
Ref. lH], performed within LIA (and after Ref. fl was published), could poten- 
tially identify the local self-crossings driven cascade. Moreover, the authors did 
observe production of vortex rings. However, they interpreted that as the pure 
Feynman's cascade. A direct numeric support for the local self-crossing scenario 
can be found in Ref. [17], where the behavior of a vortex line attached to a vibrat- 
ing object was studied within the LIA. It was observed that as a response to the 
external perturbation the vortex line gets kinky, after which frequent local self- 
crossings result in a systematic production of vortex rings. 

Still very desirable is a simulation of the zero-temperature decay of non- 
structured turbulence within the full Biot-Savart description. Such a simulation 
could quantify our very rough order-of-magnitude estimate of the crossover be- 
tween local-self-crossing-driven and pure Kelvin-wave cascades and answer the 
fundamental question of whether the inertial range of the cascade of self-crossings 
is large enough for realistic A < 15 to be taken seriously in the experimental con- 
text. 

Turning now to the Richardson-Kolmogorov regime, we emphasize the recent 
experimental progress fl9h that allowed to probe the vortex line length as a func- 
tion of temperature. As we discussed in Sec.|5l the experiment is in a very good 
qualitative agreement with our crossover scenario, the observed pronounced in- 
crease of the vortex line length at friction coefficient ^ 10~^ being indicative of 
the line fractalization. Nevertheless, it is still a circumstantial rather than direct 
evidence. A direct experimental evidence for the local self-crossings driven sce- 
nario might be an observation of small vortex rings emitted by spatially localized 
Richardson-Kolmogorov superfluid turbulence. 

On the theoretical side, a simulation of the crossover from Richardson-Kolmo- 
gorov to Kelvin-wave regime might be very instructive, especially given that nu- 
merically one can set A to be arbitrarily large to render the role of this parame- 
ter more pronounced. Yet another missing theoretical piece — in view of only an 
order-of-magnitude estimate and a technical mistake in our Ref. |[T3ll . mentioned 
in Sec. |3] — is an accurate result for the numerical prefactor in the energy flux of 
the pure Kelvin-wave cascade, Eq. ( 1431 ). 
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Summarizing our results for the crossover from Richardson-Kolmogorov to 
Kelvin-wave regime, we cannot leave without a discussion the work by L'vov, 
Nazarenko, and Rudenko lioll . which tri gge red our interest to the problem. As 
discussed in Sec. |5] the authors of Ref. 1I20I1 put forward an idea of bottleneck 
accumulation of energy in the wavenumber space between the classical-field hy- 
drodynamic modes and Kelvin waves. The notion of bottleneck in Ref. 1I20I1 is 
unequivocally defined as being associated with the absence of any channel capa- 
ble of accommodating the energy flux produced by the Richardson-Kolmogorov 
cascade, which implies thermalization of hydrodynamic modes in a certain range 
of length scales above the Kelvin-wave regime. In our Ref. iH (see also SecHl, 
we pointed out that this idea is inconsistent with the above-mentioned preemptive 
role of reconnections as an energy transport mechanism. While agreeing with us 
ll57[|22ll on the importance of the self-induced motion of vortices in the smallest — 
of the size ^ defined in Eq. (ITSl l — Richardson-Kolmogorov eddies, and attempt- 
ing to revise their theory accordingly in Ref. ll22h . the same authors still insist 
on unimportance of reconnections and support the bottleneck idea. Note, how- 
ever, that once we start from the grounds that the self-induced motion dominates 
vortex-line dynamics, reconnections are unavoidable — the lines simply do not in- 
teract with each other, and thus there is nothing to stop them from coming into 
a direct contact and reconnecting. We explicitly demonstrate in Sec. |5] that these 
reconnections serve as an efficient channel of transferring the energy flux to lower 
scales thereby ruling out a possibility of the bottleneck. 

In a certain loose sense, term 'bottleneck' might be associated with the phe- 
nomenon of fractalization of vortex lines (cf. Ref. lfl9h l since it is essentially an 
accumulation of energy in some region of the inertial range. We, however, see 
two reasons — one specific to the context of superfluid turbulence and, most im- 
portantly, a generic fundamental one — why doing so could be misleading. The 
specific reason is that the term has been already used in the context of the inter- 
face between hydrodynamic and Kelvin-wave modes, while fractalization of the 
lines is a purely Kelvin-wave effect, which happens deep (in the A ^ 00 limit) 
in the quantized regime. The fundamental reason is that the fractalization of the 
lines does not at all imply the absence of a channel(s) capable to accommodate 
the energy flux — the fractalization starting at a given wavelength scale does not 
result in any back action on larger scales, such as the development of a thermal- 
ized distribution above the bottleneck, which is a crucial ingredient of bottleneck 
scenarios. A more relevant term, putting the effect of vortex line fractalization in 
a broader context of exotic cascade phenomena, might be the deposition of the 
cascading constant of motion (e.g., energy in our case) within the cascade inertial 
range. While re-distributing the cascading quantity in a non-trivial way along the 
inertial range, the effect of deposition does not plug the cascade, thus leaving its 
generic properties intact. 
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